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ABSTRACT 

Directivity  patterns  at  3.17  MHz  and  5.1  MHz  are  calculated  for 
the  HF  antenna  array  at  the  high  power  HF  heating  facility  at  the 
Arecibo  Observatory  in  Puerto  Rico.  The  pattern  was  calculated  using 
pattern  multiplication  and  method  of  moment  techniques.  The 
calculated  pattern  is  shown  to  be  a  good  approximation  to  an 
experimentally  measured  pattern  in  one  plane  of  the  array.  A  simple 
model  was  used  to  approximate  the  effect  of  the  pattern  on  the 
frequency  response  of  ELF/VLF  signals  generated  by  the  HF  heating. 

The  frequency  response  was  determined  at  two  ELF/VLF  receiver  sites. 
Results  show  that  ELF/VLF  generated  by  side  lobes  of  the  HF  pattern 
have  sufficient  strength  to  create  a  ELF/VLF  interference  pattern  at 
receiving  locations. 


EFFECT  OF  HF  HEATING  ARRAY  DIRECTIVITY  PATTERN  ON  THE 


FREQUENCY  RESPONSE  OF  GENERATED  ELF/VLF 

INTRODUCTION 

ELF/VLF  generaCioa  experimeaCs  were  conducted  at  the  Arecibo 
Observatory  (A.O. )  in  Puerto  Rico.  The  HF  heating  faciltiy  for  A.O. 
is  located  at  18°  29'  N  and  66°  40'  W  geographic  latitude  and 
longitude  respectively.  The  ELF/VLF  receiving  site  was  located  7.7  km 
from  the  heating  facility  and  238°  to  the  east. 

The  motion  of  the  ionospheric  plasma,  in  the  presence  of  the 
earth's  magnetic  field,  causes  natural  ionospheric  currents  to  flow. 

By  changing  the  conductivity  of  a  small  portion  of  the  ionosphere,  the 
natural  currents  within  that  portion  can  be  modulated.  Modulating  the 
currents  in  Che  ELF/VLF  range  causes  a  ELF/VLF  signal  to  be  radiated 
by  the  ionosphere. 

The  conductivity  in  the  ionosphere  is  dependent  upon  the  electron 
collision  frequency,  which  is  in  turn  dependent  upon  the  electron 
temperature.  An  HF  electromagnetic  wave  is  absorbed  by  the  ionosphere. 
The  EM  wave  adds  kinetic  energy  to  the  electrons,  which  in  effect 
increases  the  electron  temperature.  Thus,  by  modulating  the  HF 
transmission  at  a  ELF/VLF  rate,  the  ionospheric  conductivities  will  be 
modulated  at  the  same  rate,  and  a  ELF/VLF  signal  will  be  radiated  from 
the  ionosphere. 

The  antenna,  radiating  Che  HF  signal  heating  the  ionosphere,  has 
a  pattern  consisting  of  a  main  beam,  side  lobes,  and  possibly  grating 
lobes.  Heating  occurs  where  each  of  these  penetrates  Che  ionosphere. 

By  determining  the  pattern  of  Che  HF  antenna,  the  spatial  distribution 
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of  Che  heacing  in  Che  ionosphere  can  be  decermined.  Thus,  Che 
locacion  in  Che  ionosphere  and  Che  inCensicy  of  each  of  che  ELF/VLF 
radiacing  scurces  can  be  decermined.  From  Chis  Che  characCerisCics  of 
Che  ELF/VLF  radiacion  from  Che  ionosphere  can  be  calculaCed. 

This  secCion  will  describe  Che  calculaCion  of  an  approximaC ion  Co 
Che  Arecibo  HF  heacing  array  direccive  gain  paCCern  and  apply  Che 
resales  Co  find  a  zero  order  approximacion  Co  an  ELF/VLF  radiaCing 
array.  The  Cechnique  employed  Co  calculaCe  Che  HF  array  paCCern  is 
one  which  uses  Che  combinaCion  of  pactern  mulciplicaCion  Cechniques 
and  compuCer  numerical  analysis.  The  numerical  program  used  was  Che 
Ancenna  Modeling  Program  (AMP)^^^.  The  AMP  ouCpuC  was  Chen  used  wich 
analytical  equations  in  a  program  written  at  the  Ionosphere  Research 
Laboratory  at  Penn  State  University  to  carry  out  the  pattern 
multiplication. 


PATTERN  MULTIPLICATION  THEORY 
The  pattern  multiplication  technique  is  based  upon  Che 
calculation  of  Che  total  pattern  of  an  array  by  caking  che  product  of 
an  array  factor  (AF )  with  the  elemental  pattern.  The  array  is  made  up 
of  identical  elements.  The  elemental  pattern  is  Che  paCtern  of  an 
individual  element  of  the  array.  The  AF  is  obtained  by  replacing  each 
of  the  elements  of  che  array  with  an  isotropic  radiator  and 
calculating  Che  pattern  for  the  array  of  isotropic  radiators.  A 
detailed  discussion  of  pattern  multiplication  can  be  found  in 
reference  (2).  A  description  of  the  theory  used  in  this  analysis 


follows. 
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For  example,  assume  thet<.  are  an  even  number  "N"  of  col  inear 
isotropic  radiators  separated  by  a  distance  "d"  as  shown  in  figure 
(1-1).  The  far  field,  at  a  point  "P",  due  to  the  n^^  radiator,  is 
proportional  to  a  complex  current  amplitude,  ,  a  phase  factor, 
e“j6Rn,  and  is  inversely  proportional  to  the  distance  from  the 
radiator  to  "P",  equation  (1-1). 

En  a  In[e“j2^n/(4TTR^)]  (1-1) 

The  far  field  approximation  states  that  "P"  is  far  enough  away 
that  "Rn"  can  be  assumed  to  be  parallel  to  "R" ,  and  the  length  of  "Rn" 
is  approximately  equal  to  "R".  Under  these  conditions  the 
approximations  in  equation  (1-2)  can  be  made. 


1/R  «  1/R(i  (l-2a) 

R^  as  (2  bl  -l)(d/2)  cos  a  +  R  ;  n  <  0  ( l-2b ) 

Rjj  a  R  -  (2  bl  -l)(d/2)  cos  a  ;  n  >  0  (l-2c) 

While  the  small  difference  in  length  of  "Rn"  can  be  neglected  in  the 
'l/Rfi"  term,  these  differences  can  have  a  significant  effect  on  the 
phase  term,  e~J6Rn.  incorporating  the  far  field  approximation 
equation  (1-2)  with  equation  (1-1),  the  total  field  at  "P"  can  be 
expressed  (1-3). 


-1  N/2 

E  "  L  En  +  E  En 

n=(H/2)  n»l 


-1 

(  E  ^n  e”j®t^2|n|-l)(d/2)cosa+R] 
n*(  -N/2) 


N/2 

+  E  W  -i^('l/2)cos  a  ]][i/(4ttR)]  (1-3) 

n^l 


Assume  that  "In"  is  equal  to  a  constant  "In"  and  collect  all  the  like 
terms.  Equation  (1-3)  reduces  to  equation  (1-4). 
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^  N/2 

E  =  (Io/47tR)  e~jSR[  ^  e~j^n  +  ^  eJ^n] 

n-(-N/2)  n=l  (1-4) 


Yn  “  8(2  kil  -l)(d/2)cos  o 

The  first  term  is  a  constant  for  a  fixed  value  of  "R".  The 
second  term  is  dependent  on  alpha.  It  defines  the  antenna  pattern  and 
is  the  AF.  The  AF  is  given  in  equation  (1-5). 

N-1  N/2 

AF  >  ^  e-JTn  +  ^  e-J^n  (1-5) 

n“(-N/2)  n=l 

Noting  that  Y-n  equal  to  Y^.  equation  (1-5)  can  be  simplified  to 
equation  (1-6). 


N/2  .  .  N/2  N/2 

AF  =  ^  (ei)^n  +  e”j)^n)  =  2^  cos  Yn  ’  cosfE  (2n-l  )(d/2)cos  a  ] 
n=l  n*l  n*l 

(1-6) 


Since  the  2  in  equation  (1-6)  is  a  constant,  it  may  be  dropped  from 
the  array  factor.  In  addition,  in  equation  (1-6)  2n-l  can  be  replaced 
with  n,  and  the  result  simplified  to  equation  (1-7). 
n-1 

AF  “  ^  cos [nB(d/2)co8  a ]  (1-7) 

n*l ,3,5. . . 

Equation  (1-7)  is  the  array  factor  for  a  col inear  array  of  an  even 
number  of  N  isotropic  radiators  with  equal  current  amplitudes.  The 
angle  a  is  measured  from  the  line  of  the  array  to  the  observation 
point. 
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ANTENNA  MODELING  PROGRAM  (AMP)  THEORY 


The  Antenna  Modeling  Program  (AMP)  was  used  to  analyze  the  array 
element.  This  program  was  developed  by  MBA/ In format  ion  Systems. 

The  computer  program  applies  the  method  of  moments  to  the  thin  wire 
approximation  of  the  integral  equation  for  the  electric  field  due  to  a 
volume  current  distribution,  equation  (1-8).^^^ 

E(F(j)  ■  (0  J(T)*  (G(F,FQ)dv  (1-8) 


G(F,Fo)  =  -(l/4iT)[I+(l/k2)VV]g 
g  -  (e-il*  /  F-Fol  ) 
k  “  wi/uoeQ  ,  I  *  unit  2”*^  rank  tensor 

(F-FqI  “  distance  measured  from  wire  axis  (source  point) 
to  observation  point  on  the  surface. 

The  thin  wire  approximation  requires  that  the  diameter  of  the 

wire  be  small  compared  with  the  wavelength.  Thus  azimuthel  current 

flow  around  the  wire  can  be  neglected  and  the  volume  integral  in 

equation  (1-8)  can  be  changed  to  a  line  integral,  equation  (1-9).^^^ 


-Sg'E^CFo)  »  (-i(i)Uo/4Tr)  ^I(s)[3*SQ-(l/k2)  (32/9s3S(j]g(r ,ro)ds 


(1-9) 


3  >  unit  tangent  at  source  point 

Sg  >  unit  tangent  at  observation  point 
I  *  (■ira2j)2iTa 


a  *  wire  radius 
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Included  in  equation  (1-9)  is  also  the  boundary  condition  for  a  metal 
surface,  equation  (1-10). 

E^can  +  E^tan  *  0 

Elcan  =  Tangential  component  of  incident  electric  field 

E^tan  “  Tangential  component  of  scattered  electric  field 

AMP  solves  equation  (1-9)  numerically  by  converting  it  into 
matrix  form.  This  is  accomplished  by  expanding  the  unknown  currents, 
I,  in  terms  of  a  set  of  basis  functions,  I,),  and  taking  the  inner 
product  of  both  aides  of  equation  (1-9)  with  a  set  of  weighting 
functions  w^,.  A  general  discussion  on  this  method  of  solution  can  be 
found  in  reference  (3). 

Equation  (1-11)  is  obtained  by  expressing  equation  (1-9)  in 
operational  format,  where  the  operator  L^p  denotes  the  integral  and 
"<A,B>"  denotes  the  inner  product  of  quantities  A  and  B. 

N 

^  An  <Wjj,  Inp  Ifi^  "  e1>  (1-11) 

n  N 

where  I  -  A^ 
n 

Equation  (1-11)  must  be  true  for  each  and  thus  may  be  written 
in  matrix  form  as  expressed  in  equation  (1-12). 


^^1  >  Iqjjp 

^'^2*  "^^2*  ^>p  ^2^  l^p 

fi 

ii  fi  II 

ti  fi  ft 

If  fi  fi 

Al 

A2 

If 

II 

If 

II 

= 

<wj ,  e1> 

Cw2)  e1> 

II 

II 

It 

ft 

^5p  ^p  ^2^ - Lop 

<W{} ,  e1^ 

(1-12) 
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Since  E^,  Ijj  and  Wj,  are  known,  by  matrix  inversion  the  values  of 
A{,  can  be  calculated. 

Specifically,  AMP  uses  sine  and  cosine  functions  as  basis 
functions  and  employs  a  method  of  collocation,  or  point  matching,  by 
choosing  the  weighting  functions  as  6  functions. 

THEORY  APPLIED  TO  A.O.  ARRAY 

To  apply  the  techniques  of  AMP  and  pattern  multiplication  to  the 
A.O.  heating  array,  the  array's  physical  characteristics  must  be 
known.*  The  HF  antenna  consists  of  a  4x8  array  of  radiating  elements. 
This  array  is  oriented  as  shown  in  figure  (1-2).  Each  of  the  elements 
in  the  array  is  constructed  in  the  shape  of  an  inverted  pyramid  with 
four  sides.  The  faces  of  Che  pyramid  are  at  an  angle  of  45°  with  Che 
ground  and  contain  two  nonplanar  log-periodic  antennas  (NLPA).  One 
NLPA  is  contained  in  Che  north  and  south  faces,  and  the  other  is 
contained  in  Che  east  and  west  faces.  A  sketch  of  an  array  element  is 
shown  in  figure  (1-3).  Note  that  the  elements  in  Che  south  and  west 
faces  have  been  rotated  180°  about  the  corresponding  face's  NLPA 
elements  feed  lines.  A  diagram  looking  down  at  Che  top  elements  of 
the  pyramid  is  shown  in  figure  (1-4). 

Both  NLPAS  are  designed  with  a  t  of  .88.  The  dimensions  of  the 
array  elements  in  the  north  and  south  faces  are  shown  in  figure  (1-3). 
The  dimensions  of  Che  east  and  west  faces  are  scaled  to  of  Che 

north  and  south  faces.  This  will  result  in  right  hand  circular 
polarization  radiation  when  the  north  and  south  faces  are  fed  180°  out 

♦Note:  See  appendix  III  for  additional  information  on  HF  antenna 

geometry. 
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Figure  1-2  HF  heating  array 


Figure  1-3  HF  heating  array  element 


Figure  1-A  View  of  top  elements  looking  down  at  pyramid 


ALL  WIRE  RADII=2  0  mm 
ALL  DIMENSIONS  IN  METERS 
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Figure  1-5  Non-planar  log-periodic  antenna  semi-structure  dimensions 
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of  phase  with  the  east  and  west  faces.  When  fed  in  phase,  left  hand 
circular  polarization  radiation  will  be  transmitted. ^  The  north  and 
south  and  the  east  and  west  faces  are  fed  against  ground  by  separate 
transmitters.  This  gives  the  capability  to  adjust  the  phase  between 
different  faces  and  thus  change  the  radiation  polarization. 

The  array  element  was  analyzed  using  the  AMP  computer  program. 
Then  the  pattern  multiplication  technique  was  used  to  calculate  the 
total  array  pattern. 

The  array  element  was  analyzed  using  AMP  for  frequencies  of 
3.17  MHz  and  5.1  MHz.  The  data  file  used  is  listed  in  Appendix  I 
program  1.  The  GW  and  GM  cards  generate  the  antenna  structure 
depicted  in  figure  (1-3)  and  orient  it  with  respect  to  the  coordinate 
axis  shown  in  figure  (1-4).  The  GN  card  specified  the  conductivity 
(.03O/m)  and  relative  permativity  (20)  of  the  ground  below  the  HF 
array.  The  actual  conductivity  and  permativity  for  the  heater  site 
was  not  known,  but  based  on  the  fact  that  maps  indicate  the  heater 
array  is  located  on  marshy  ground,  a  conductivity  and  a  relative 
permativity  for  "good"  earth^^^  were  used.  As  shown  in  the  RP  card, 
the  power  gain  was  computed  for  2.5  degree  steps  in  "theta"  and  5 
degree  steps  in  "phi." 

In  order  to  use  the  output  of  AMP  in  the  pattern  multiplication, 
it  was  necessary  to  develop  an  elemental  pattern  function.  Given  a 
"theta"  and  "phi",  the  function  returns  a  value  of  the  power  gain  in|» 
that  direction.  To  accomplish  this,  the  power  gains  form  the  AMP 
output  for  selected  constant  "phi"  surfaces  were  combined  with  a 
linear  interpolation  scheme.  The  selected  values  of  "phi"  are  listed 


in  table  I. 
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Figures  (1-6)  and  (1-7)  show  the  comparison  of  the  interpolated 

values  (shown  by  X)  with  the  AMP  results  (shovai  by  •).  The  maximum 

error  is  approximately  one  db.  Figure  (1-6)  is  a  plot  of  power  gain 

as  a  function  of  "phi"  for  a  constant  "theta."  The  power  gain 

decreases  as  Che  radius  of  the  polar  plot  increases.  The 

unsytnmetrical  nature  of  the  plot  is  due  to  the  unsymmetr ical  nature  of 

the  array  element.  Figure  (1-7)  is  a  plot  of  power  gain  as  a  function 

of  "theta"  for  a  constant  "phi." 

3.17  MHz  5.1  MHz 

Phi  (deg)  Phi  (deg) 


0 

10 

90 

60 

130 

90 

180 

no 

210 

150 

270 

180 

320 

215 

360 

270 

320 

Table  I.  Selected  Values  of  Phi  for  Interpolation  Routines 

The  total  array  pattern  for  the  A.O.  array  shown  in  figure  (1-1) 
was  calculated  by  taking  Che  product  of  two  separate  array  factors. 

In  essence  this  is  a  pattern  multiplication.  The  array  pattern  of  one 
array  becomes  the  elemental  pattern  of  the  ocher  array.  One  array  is 
the  4-element  array  on  the  north  and  south  line.  The  other  array  is 
an  8-e lament  array  on  an  east  and  west  line. 

The  expansion  of  equation  (1-7)  for  the  4-element  and  8-element 
arrays  can  be  simplified  by  using  trigonometric  identities.  These 
trigonometric  identities  were  obtained  from  Chebysheff  polynomials,  as 
shown  in  Appendix  II.  The  simplified  equation  for  Che  antenna  factors 
are  given  in  equations  (1-13)  and  (1-14).  The  4-element  array  factor 
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Figure  l-6a  Power  gain  vs.  phi  for  a  constant  theta,  3.17  MHz 
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Figure  l-6b  Power  ?ain  vs.  ohi  for  a  constant  theta,  5.1  MHz 
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Figure  l-7a.2  Power  gain  vs.  theta  for  phi=40 
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Figure  l-7a.5  Power  gain  vs.  theta  for  phl=280 


Figure  l-7b.l  Power  gain  vs.  Clieta  for  phl=0  ,  5.1  Mllz 
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Figure  l-7b.3  Power  gain  vs.  theta  for  phi=250 
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is  equation  (1-13),  and  the  8-element  array  factor  is  equation  (1-14). 
Figure  (1-8)  shows  the  orientation  of  the  4-  and  8-eleraent  arrays  on 
an  X,Y,2  coordiate  system.  The  Y-axis  was  taken  to  be  north  and  the 
X-axis  as  east. 

AF^  *  [sin(48(d/2)  cos Cj I/3in(8 (d/2)co8  C^)  (1-13) 

AF2  =  [ain(86(d/2)cosC2)l/sin(8(d/2)co8  C2)  (1-14) 

The  array  pattern  of  the  8-element  array  is  solid  of  revolution 
about  the  X-axis,  and  the  4-eleinent  array  pattern  is  a  solid  of 
revolution  about  the  Y-axis. 

Using  the  transformation  in  equation  (1-15),  AF^  and  AF2  can  be 
transformed  to  spherical  coordinates.  The  total  array  factor  AF , 
equation  (1-16),  is  the  product  of  AFj  and  AF2. 

cos  *  sinQsin  <p  (1-15) 

cos  ^2  *  3  in  0  cos 

AF  *  AF  X  AF  =  9int^6(d/2)3ia03in<))]  3in[8S(d/2)3in0co3$]  (1-16) 

^  2  sin  [g  (d/2)sin0sin(()]  sin[B(d/2)  sin0coS(j] 

The  total  array  power  pattern  can  be  calculated  by  taking  the 
square  of  the  total  array  factor,  AF,  and  multiplying  it  by  the 
elemental  pattern  function,  which  is  an  interpolation  of  the  AMP 
results.  To  achieve  the  goal  of  the  calculation  of  a  directive  gain 
pattern,  a  correction  factor  must  be  determined.  This  factor  is  a 
result  of  neglecting  the  constants  in  calculating  the  array  factor. 

The  directive  gain  "in  a  given  direction  is  defined  as  the  ratio 
of  the  radiation  intensity  in  chat  direction  to  the  average  radiated 
power. Since  the  constants  which  were  neglected  in  calculating  AF 


are  also  contained  in  Che  calculation  of  the  average  radiated  power 
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and  the  directive  gain  is  a  ratio,  these  same  constants  must  also  be 
neglected  in  calculation  of  the  average  radiated  power.  Thus  the 
calculation  of  the  average  radiated  power  becomes  the  correction 
factor  necessary  to  convert  the  array  power  pattern  to  an  array 
directivity  pattern. 

Equation  (1-17)  was  used  to  calculate  the  average  radiated  power, 


-  ■ 


2ir^ir/2 

\  [AF(0,<j>)]2  X  10  (ELF(0,(ji)/lO)  3in0d6  dd 

0  Jo  4it 


(1-17) 


ELF(0,  4i)  =  Elemental  Power  gain  from  interpolated  AMP  output 

AF(  0,  (p)  =  Total  Array  Factor 

The  integration  was  performed  numerically  using  Simpson  integration, 
equation^  (1-18).  The  programs  are  given  in  Appendix  I,  programs  2 
and  3. 

^  f(x)dx  =  (Ax/3)  [f(xo)+4f(xi )+2f(x2)+4f(x3)+2f(x4)+. . .  (1-18) 

+2f (x2n-2  (*2n-l )+f(x2n) ^ 

Ax  “  (b-a)/(2n) 

The  first  quardrant  integration  was  carried  out  for  two  cases. 

One  case  was  with  a  "phi"  step  size  of  2.5  degrees;  the  other  case  was 
with  a  "phi"  step  size  of  1  degree.  A  "theta"  step  size  of  1  degree 
was  used  in  both  cases.  No  significant  difference  was  found  in  the 
result  of  the  integrations.  Based  on  this  result,  the  step  sizes 
chosen  for  the  total  integral  were  1  degree  and  2.5  degrees  for 
"theta"  and  "phi"  respectively.  The  correction  factors  determined 
were  9.82  db  and  8.62  db  for  5.1  MHz  and  3.17  MHz  respectively. 
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Equation  (1-19)  combines  the  total  array  factor  AF,  the  elemental 
pattern  ELF,  and  the  correction  factor  to  calculate  the  directive  gain 
pattern  for  the  A.O.  heating  array. 

D(0,i^)  =  20  log[AF(  0,  ij))  ]  +  ELF(0,(j))  -  Correction  factor  (1-19) 

Figure  (1-9)  is  a  plot  of  the  pattern  in  the  "phi"  equal  zero 
plane  (north-south  plane).  The  "x's"  are  experimentally  measured 
values.^  The  values  were  measured  from  a  Boeing  707  aircraft  at  2900 
ft  (3.84  km).  The  plane  was  flown  over  the  A.O.  array  on  a 
north-south  line,  while  the  heater  was  operating  at  5.1  MHz.  The  plot 
shows  that  the  array  pattern  obtained  by  the  combination  of  pattern 
multiplication  and  numerical  techniques  is  a  good  approximation  of  the 
A.O.  array  pattern. 

Figures  (1-10)  and  (1-11)  are  plots  of  the  directivity  pattern  of 
the  A.O.  heater  array  for  3.17  MHz  and  5.1  MHz  respectively.  Each 
plot  is  the  variation  of  the  directive  gain  with  "theta"  in  a  constant 
"phi"  plane.  "Phi"  is  varied  in  5  degree  steps  from  0  degrees  to  180 
degrees.  Two  additional  "phi"  plane  patterns  have  been  plotted  in 
figures  (1-12)  and  (1-13).  These  two  figures  are  the  directivity 
patterns  in  the  "phi"  equal  121.5*^  and  146°  planes  respectively. 

These  patterns  are  in  planes  corresponding  to  the  direction  of  Los 
Canos  and  the  A.O. ,  respectively.  The  programs  used  to  make  the 
directivity  patterns  are  given  in  Appendix  I,  programs  4  and  5  for 
3.17  MHz  and  5.1  MHz  respectively. 
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Figure  1-11  Directive  gain  pattern  for  Arecibo  Observatory 
HF  heating  array.  Frequency  “  5.1  MHz. 
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Gracing  lobes  (lobes  which  have  Che  same  incensicy  as  Che  main 
beam^)  in  Che  anCenna  facCor  will  occur  when  boCh  Che  numeraCor  and 
denominaCor  of  boCh  Cerms  of  equaCion  (1-15)  are  zero.  This  will 
occur  when  equaCion  (1-20)  is  saCisfied. 

g(d/2)  sinQsiniJ)®  0  or  it 

(1-20) 

6(d/2)  sinQcosiJi  =  0  or  ir 

For  Che  A.O.  array  3d/2  is  equal  Co  2.82  and  4.54  for  3.17  and 
5.1  MHz  respecCively.  Since  sin0  sini^  and  sinG  cos())  are  never  larger 
chan  1,  Che  3.17  MHz  paCCern  does  noC  and  should  noC  have  gracing 
lobes.  However,  gracing  lobes  will  be  presenC  in  Che  5.1  MHz  paCCern 
because  4.54  is  larger  chan  rr. 

These  gracing  lobes  will  occur  for  Che  angles  given  in  Cable  II. 
As  can  be  seen  in  Che  plocs  of  figure  (1-11)  major  lobes  do  occur  aC 
Chese  angles.  They  are  aCCenuaCed  when  Chey  are  mulciplied  by  Che 
elemenCal  paCCern  during  Che  calculacions  of  Che  CoCal  array  pacCern. 

9(deg)  (<>(deg) 


43.8 

0 

43,8 

90 

43.8 

180 

43.8 

270 

78.1 

45 

78.1 

135 

78.1 

225 

78.1 

315 

Table  II.  LocaCion  of  Gracing  Lobes  in  5.1  MHz  paCCern 
ELF/VLF  ARRAY  MODEL 

Having  esCablished  a  direcCive  gain  paCCern,  ic  remains  Co  relaCe 
Che  paCCern  Co  Che  heaCing  of  Che  ionosphere.  Richardson®  has  shown 
chac  Che  largesc  change  in  conduce ivicy  caused  by  che  heacing  occurs 
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at  approximately  a  70  km  altitude.  As  a  zero  order  approximation,  the 
heating  pattern  can  be  projected  on  a  plane  located  at  a  70  km 
altitude.  The  location  and  relative  intensity  of  the  major  heated 
regions  can  be  found.  By  placing  elementary  dipoles  with  the  same 
relative  amplitude  of  current  at  the  respective  heated  regions,  a 
field  intensity  at  a  receiving  site  on  the  ground  can  be  calculated. 

A  correction  factor  is  needed  to  project  the  pattern  onto  a  70  km 
altitude  plane.  The  pattern  shows  the  relative  distribution  of  the 
power  on  a  sphercial  surface  of  radius  "R. "  Since  the  distance  to  a 
plane  increases  when  "theta"  is  greater  than  zero,  the  power  density 
on  the  plane  will  decrease  from  that  indicated  by  the  pattern.  The 
power  is  being  spread  over  a  larger  spherical  surface  as  "R"  is 
increased.  The  projection  of  the  pattern  on  a  plane  surface  requires 
multiplying  the  pattern  by  an  attenuation  factor  of  cos^©.  This 
attenuation  is  plotted  as  a  function  of  "theta"  in  figure  (1-14). 

Programs  6  and  7  in  Appendix  I  were  used  to  compute  a  pattern  on 
a  square  section  of  a  plane  (122  km  north  and  south  by  122  km  east  and 
west  of  the  main  beam)  at  a  70  Ion  altitude.  The  data  were  then 
plotted  using  Statistical  Analysis  System  (SAS).^^>^®^  Figures  (1-15) 
and  (1-16)  show  the  relative  power  density  on  a  70  Ion  altitude  plane 
for  3.17  MHz  and  5.1  MHz  respectively.  Only  levels  above  that  of  an 
isotropic  radiator  (0  db)  were  plotted. 

The  122  km  dimension  of  the  plotted  data  corresponds  to  an  angle 
of  "theta"  equal  to  60.2  degrees.  The  major  lobes  created  by  the 
antenna  factor  grating  lobes  at  78.1  degrees  are  not  seen  in  figure 
(1-16).  These  lobes  are  for  the  most  part  attenuated  to  the  zero 
reference  level  by  the  long  propagation  path.  Table  III  gives  the 
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Fii?ure  1-14  Propagation  loss  due  to  power  spreading 
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location  and  relative  power  density  of  these  lobes  on  the  70  km 
altitude  plane. 

Data  on  the  lobes  shown  in  figures  (1-15)  and  (1-16)  are  given  in 
tables  IV  and  V.  The  current  was  assumed  to  be  proportional  to  the 
square  root  of  the  power  density.  To  find  the  average  current  per 
unit  length  for  each  lobe,  the  current  was  integrated  over  the  area  of 
the  plane  disturbed  by  the  lobe  .»nd  divided  by  the  length  of  the 
region.  The  current  distribution  on  the  plane  was  approximated  using 
a  pyramid  with  a  quadrilateral,  as  a  base.  The  length  and  width  of  the 
heated  region  are  the  diagonals  of  the  quadrilateral  base;  the  peak 
current  is  the  altitude.  The  volume  of  the  pyramid  is  equal  to  one 
third  of  the  area  of  the  base  times  the  altitude.^®  The  area  of  the 
quadrilateral  base  is  one  half  of  the  product  of  the  diagonals  times 
the  sine  of  the  angle  between  the  diagonals. So  the  formula  for 
calculating  the  average  current  per  unit  length  becomes  equation 
(1-21). 

lav  =  (l/3)[(l/2)a-b  sin  90l ( Ip/a )=( 1/6 )  bip  (1-21) 

a  *  length  (diagonal) 
b  *  width  (diagonal) 

Ip  =  peak  current  (altitude) 

In  this  zero  order  approximation  each  of  the  lobes  is  represented 
as  an  elementary  current  element  source.  The  source  has  a  current 
equal  to  the  average  current  calculated  in  equation  (1-21),  a  length 
equal  to  the  length  of  the  lobe,  and  is  located  at  the  x,y  coordinates 
of  the  peak  power  density  for  the  lobe  cn  the  70  km  altitude  plane. 

The  location,  length,  and  average  current  are  summarized  in  tables  IV 


and  V. 
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Pd 


Locat ion 
x(kia)  y(km) 

Re  lac ive 
Power 
Dens iCy 
(db) 

a 

Length 

(km) 

b 

Width 

(km) 

Ip 

Peak 

Current 

(lOpd/20) 

lav 

Average 
Current 
(1/6  a*b‘Ip)a 

-63 

0 

3.022 

24 

12 

1.416 

2.832 

-30 

0 

10.87 

22 

18 

3.495 

7.485 

0 

0 

24.983 

40 

20 

17.748 

59.16 

30 

0 

10.215 

22 

16 

3.242 

8.645 

62 

0 

1.295 

15 

7 

1.161 

1.355 

0 

-38 

1.624 

19 

7 

1.207 

1.408 

0 

-25 

5.875 

26 

9 

1.967 

2.951 

0 

-14 

11.239 

32 

9 

3.647 

5.471 

0 

14 

11.239 

32 

9 

3.647 

5.471 

0 

25 

5.875 

26 

9 

1.967 

2.951 

0 

38 

1.624 

16 

7 

1.206 

1.407 

Table 

IV.  Lobe 

SCaCistica 

for  Frequency  “  3.17 

MHz 
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Locat ion 

X  ( km )  y  ( km ) 

Re lat ive 
Power 

Dens ity 
(db) 

a 

Length 

(km) 

b 

Width 

(km) 

'P 

Peak 

Current 

(lOpd/20) 

lav 

Average 

Current 

(1/6  a‘b'Ip)/a 

-60,5 

0 

15.519 

54 

14 

5.970 

13.93 

-34.0 

0 

10.141 

14 

14 

3.214 

7.499 

-18.0 

0 

12.205 

12 

10 

4.076 

6.793 

0 

0 

24.113 

24 

12 

16.056 

32.113 

18.0 

0 

11.134 

12 

10 

3.603 

6.006 

33.0 

0 

8.  306 

14 

10 

2.530 

4.217 

59.0 

0 

11.279 

43 

14 

3.664 

8.549 

0 

-65 

14. 301 

26 

28 

5.189 

24.213 

0 

-48 

5.331 

15 

8 

1.847 

2,463 

0 

-38 

3.715 

12 

6 

1.5  34 

1.5  34 

0 

-30 

3.586 

12 

5 

1.511 

1.259 

0 

-22 

4.825 

14 

5 

1.743 

1.452 

0 

-15 

6.969 

16 

5 

2.231 

1.859 

0 

-9 

11.079 

20 

6 

3.581 

3.581 

0 

9 

11.079 

20 

6 

3.581 

3.581 

0 

15 

6.969 

17 

5 

2.231 

1.859 

0 

22 

4.825 

15 

5 

1.743 

1.452 

0 

30 

3.586 

14 

5 

1.511 

1.259 

0 

38 

3.715 

14 

6 

1.534 

1.534 

0 

48 

5.331 

18 

8 

1.847 

2.463 

0 

65 

14.301 

28 

28 

5.189 

24.213 

Table  V 


Lobe  Statistics  for  Frequency  =5.1  MHz 
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The  magnetic  field  due  to  a  current  element  is  given  in  equation^ 
(1-22).  The  source  is  located  at  the  origin  and  along  the  Z'  axis. 

is  the  magnetic  field  intensity  in  the  source  coordinate  system 
X'  ,Y'  ,Z'.  Since  is  a  vedtor  and  there  are  a  number  of  sources  at 

different  locations  whose  fields  need  to  be  superimposed  at  the 
observation  point,  it  would  be  beneficial  to  translate  the  fields  to 
the  observation  frame  of  reference,  X,Y,Z. 

-  ( (IdD/Amr' ]ain0'e~j^r' [j6  +  (1/r')]  (1-22) 

The  source  and  observation  coordinate  systems  are  shown  in  figure 
(1-17).  The  source  is  located  at  point  O(X0,Y0, ZO)  and  along  the 
Z'-axis,  which  is  parallel  to  the  X-axis  in  the  observation  point 
coordinate  system.  The  Y-  and  Y'-axis  are  also  parallel  and  the 
X'-axis  is  in  the  negative  z  direction.  From  the  geometry  of  the 
figure  the  relationships  given  in  equation  (l-2j)  can  be  found.  In 
addition,  using  the  transformation  from  spherical  coordinates  to 
Cartesian  coordiantes,  equation^2  (1-24),  the  expression  for  the 
magnetic  field,  equation  (1-22),  can  be  transformed  to  the  observation 
point  coordinates. 


X'  - 

XP-XO 

Y*  - 

YP-YO 

Z*  “ 

ZP-ZO 

Y0)2+( ZP-ZO )2) 1/2 

r'  - 

(XP-XO )2+(yP- 

sin©' 

-  ( (ZP-ZO ))2 

+( YP-YO )2] 1/2 

Hy  »  Hy'  -  (x' )/Vx'^+y'^ 

-Hg  -  Hx'  -  (y'  )/Vy'2+x'2  (1-24) 

Hx  -  Hx'  -  0 

The  phase  term  e~j6*'  is  important  when  performing  the 
superposition  of  the  fields  from  all  of  the  sources  at  the  observation 
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OBSERVATION  COORDINATES  X,Y,Z 
SOURCE  COORDINATES  X',Y‘,Z’ 


Figure  1-17  Relative  orientation  of  observation  and  source  conrdina 


I  /2 
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point.  Small  differences  in  Che  value  of  r'  can  be  magnified  by 
causing  a  significant  change  in  Che  term's  value.  Thus  it  would 
affect  Che  summation  of  Che  real  and  imaginary  parts  of  all  the 
sources. 

An  additional  term  must  be  added  Co  8r'  Co  maintain  Che  proper 
relationship  between  Che  sources.  This  is  because  Che  HF  heating 
pulse  must  travel  different  path  lengths  to  the  source  region.  The 
difference  in  path  lengths  causes  a  time  delay  between  the  sources, 
equation  (1-25).  This  delay  can  be  expressed  in  terms  of  an  additonal 
path  length  AR,  equation  (1-26).  This  allous  Che  phase  term,  e'  -jBr' ^ 
in  equation  (1-22)  Co  be  expressed  as  +AR)^  By  combining  the 

phase  delay,  equations  (1-23)  and  (1-2A)  inro  equation  (1-22),  an 
expression  for  Che  magnetic  field  intensity  at  the  point  of 
observation  and  in  terms  of  Che  observation  point  coordinate  system 
can  be  found,  equation  (1-27). 


Time  delay  =*  ( ZO^-KYO^-t-XO^ )  1/2-70 

C 


(1-25) 


Phase  delay  -  (ZO^-i-YO^-t-XO^  )l/2-70  -  (Ztt/X)  [  (Z02+Y02+X02  )l/2-70] 

C 

-  BAR  (1-26) 


if  >  Idl[(ZP-Z0)2-*-(YP-Y02)l/2 

4nl(  XP-XO  )  •^+ ( YP- YO )  •^ + (  ZP-ZO  ^ 


.-jB[(XP-XO)2+(YP-YO)2  + 


(ZP-ZO)^)^^^  +  (ZP^  +Y0^+X0^)^'^^-70 


(YP-YO) 

((XP-X0)2+(YP-Y0)^)J 


r  (xp 

|{XP-X0)2+ 


-XO ) 

( YP-YO  )'^)J 


j8+  _  _  _ 1  _ _ 

[ (XP-XO )^+(YP-Y0 )2+( ZP-ZO )2 ) ] J 


(1-27) 
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Programs  8  and  9  Appendix  I  were  written  to  calculate  the 
strength  of  the  magnetic  field  at  an  observation  point  due  to  the  two 
source  patterns  of  figures  (1-15)  and  (1-16).  For  each  of  the  two 
cases  a  current  element  was  placed  at  the  location  of  the  lobes,  and  a 
total  magnetic  field  was  calculated  at  an  observation  point 
corresponding  to  Los  Canos.  The  calculations  were  carried  out  for  the 
frequencies  listed  in  table  VI.  These  frequencies  correspond  to  the 
frequencies  used  during  the  ELF/VLF  experiments  conducted  in  Puerto 
Rico.  Figure  (1-18)  shows  the  result  of  the  calculation.  It  is  a 
plot  of  relative  magnetic  field  strength  as  a  function  of  frequency. 


Frequency 
(kHz  ) 


Frequency 
(kHz ) 


.479386 

1.0 

1.25 

1.506024 

2.0 

2.293578 


2.5 

2.793296 

3.144654 

3.448276 

4.0 

4.464286 

5.0 


Table  VI.  Experimental  Frequencies 


In  order  to  determine  the  effect  of  the  lobes  on  the  value  of  the 
magnetic  field  at  the  observation  point,  a  calculation  was  done  with 
only  the  main  lobe  acting  as  the  source.  The  plot  of  the  relative 
field  strength  as  a  function  of  frequency  is  given  in  figure  (1-19). 

A  comparison  of  figures  (1-18)  and  (1-19)  shows  that  the  lobes  have  a 
significant  effect  on  determining  the  frequency  response  of  the 
ELF/VLF  radiating  source.  The  location  of  the  relative  maximums  and 
minimums  in  both  5.1  MHz  and  3.17  MHz  generated  ELF/VLF  response  have 
shifted  in  frequency.  A  significant  reduction  in  the  field  strength 


occurs  between  2.5  to  4.5  kHz  for  the  5.1  MHz  generated  Y  component. 


RELATJVE.  H  FIELD  STRENGTHS 
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Figure  1-19  Main  beam  current  element  VLF/ELF 
response 
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In  addition  to  computing  the  relative  field  strengths  for  the 
test  frequencies,  a  calculation  was  made  in  100  Hz  frequency  steps 
from  from  500  Hz  to  5  kHz.  The  result  is  plotted  in  figure  Cl-20). 

As  can  readily  be  seen  from  the  figure,  two  deep  minimums  occur  for 
the  5.1  MHz  generated  ELF/VLF.  One  occurs  between  600  Hz  and  700  Hz, 
and  the  other  occurs  between  3400  Hz  and  3600  Hz.  The  3.17  MHz 
generated  ELF/VLF  has  only  one  deep  minimum,  which  is  located  between 
700  Hz  and  900  Hz. 

A  calculation  was  made  for  an  alternate  receiving  site.  This 
site  is  located  at  Salinas,  Puerto  Rico,  17.98°  N  and  66.30°  W 
geographic  latitude  and  longitude  respectively.  The  frequency 
response  results  for  the  two  HF  heating  patterns  are  shown  in  figure 
(1-21).  Comparison  of  figures  (1-20)  and  (1-21)  shows  the  changes  in 
the  response  due  to  the  different  propagation  paths.  The  first 
striking  difference  is  the  disappearance  of  the  null  in  between  3  and 
4  KHz  in  the  5.1  MHz  pattern  generated  VLF.  Second,  the  relative 
maximums  and  minimums  for  the  Salinas  location  have  been  shifted  to  a 
lower  frequency.  In  general,  the  relative  field  strengths  are  lower 
for  the  Salinas  site  due  to  the  increase  in  propagation  distance. 

CONCLUSION 

This  completes  the  description  of  the  zero  order  approximation. 

In  summary,  a  pattern  for  the  Arecibo  Observatory  has  been  calculated 
using  the  technique  of  pattern  multiplication  and  AMP.  The  calculated 
pattern  for  the  "PHI"  equal  "0"  plane  compares  with  the  experimentally 
measured  pattern.  The  pattern  shows  that  there  is  enough  power  in  the 
side  lobes  and  grating  lobes  to  cause  significant  heating  at  a  70  km 
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RELATIVE  FIELD  STRENGTH 


I Igure  1-21  Current  element  array  Vl.K/ELF  response  at  Salinas.  V  component  of  magnetic  field 
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altitude.  Using  the  HF  antenna  pattern,  a  zero  order  approximation  of 
the  model  for  the  ELF/VLF  radiating  system  was  determined.  The 
frequency  responses  of  the  ELF/VLF  radiation  system  for  observation 
points  corresponding  to  Los  Canos  and  Salinas  were  calculated.  The 
calculations  determined  that  the  regions  heated  by  the  side  lobes  and 
grating  lobes  have  a  significant  effect  on  the  strength  of  the 
received  signal. 
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,  ^  .  1^  •  .  W  ,  s 

1 ,  *  0 , ,  ,  ‘1 ,  /,  F ,  .  A  7 

»  S  f  *» (1  ^  *  *1 ,  s  * 

•  •' 

,f.n, 

» » ^  t 

Ml.  t  /  f'  . .  f'l 

1 .  .  7  7  .*)  .Mo  , .  M.i  , 

.  t  »W  ,‘H»  * 

llA  r  A  (I  /  7  B  A  .M  1  .  1 

.  ,  i,  / 

.  t  /  A 

. MS ,  « , MW 

,  1  .•••»,  1 , r.  / , 

1  ,  M  ,  1  ,  f.  ,  1  , 

.  I."!'  .  i.  i.- ,  '  .  ' 

'•  1  f 

.  t  ' 

,  • 

II  ,  ‘..Ml,  1 

,  W  •  1  •  M"  iM 

1  '  .  'M  ,  <  ,  1  , 

1  A.  P*  .  A  .  AA  .  A  .  7  A  .  A 

.mm,  1 

,  t 

•  S  ,  I  .  M  , 

1  ,  '1  1  f.  , 

1  .  A'l  « .'  •  7  Ano  , ,‘i  \  . 

1  •*•/.  •  4  .  \  . 

M  .# 

II  , 

f>  , ,1.11  ,  (.mi 

*  '  .f'  ,  '  .  /w  ,  ■ 

1  ,  A7.,  1  ,  lA,  1 

r*  ♦  w  .’  t  * 

I  , 

fl,  1 

1 .  '1  f ,  1 

,  n  1 ,  1 ,  ,  1 , 

IS  TO  ,  1  ,')  7  ,  1  .IT'.  ,  1  , 

<  1  •  1  •  /n  , ,  M  t 

*  1  . 

1  , 

If.  ,  1  ,  , 

.  1 . • 

■ 

1  .  1 . 1  ,  .  "I  .  . 

♦ 

.  • 

1  ,  ,  •  , 

,  / 1 ,  .  '  , ,  ’ 

'  1  »  ,  ■ "  .  ,  1  • 

1  •  .  7f  A  ,  -  .  1  A  .  -  1  ,  A  1 

» -  f ' .  1 

It 

.  • .  .  1 

.  •  * .  .  1  •» . 

1  .  <  ' . 
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i  'STiis  • '  ’ "  r  r.'’ ^ .V7 .  - 1  s .  AA .  -  I  . . . 

_L,; ‘ ‘ A I .-.  w 

flATA  - iv  k,  -J 

ii»i 
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— rrnn — FtrAcr-TT^FT  Ji^ATrrmaTTT  T7  ta  top  t  a7  j  rr-Ainr  rATJrn 

IH  Ml  .t-U.  Ml  UK  I  I  F  ( h  ,H0II  I  i-ni 
HW  ac.Ks  (  km  1  -  Ai'lll  M  I  I  I  ■'  t  akhI  I  I  ■<  I  -  akhI  I  I  1  I  1 
IFM  <  1(1 1  I  4=1  I 

FF)=HKArl'»(FF(l4,.nl-FF(M.Jlll*F'-MI..<M 

_ fcFgsFKArKAIFK  M,.I4I-FF(  I1.J4II-4FKII.JA) _ ^ _ 

FLF=t^KACT»(FF?-FFt  l+FFt 
FLHF.=Pt>/U>. 

_  FLFi  I  II.  KSFI.FF _ _  _  _  _ _ 

r.  wS  ItF  r  Ton  I  FLF  .A~|  HfT'a  (  j  4  i  ,  A  i  HF  I  A  i  j  fi  .  AKM  I  (  I  4  I.  AKHI  (  fl  )' 

r.  7n(»  fmkmati'  ' .  ■  fi.fmfmt  ai,  facIm»=  •  ,fa.?.  •  mh  • ,  a|  ax  .fa  .  ?  )  l 
HOO  ’  FDKMAK'  '.Fin, a,'  mhi  Mill  IN  kan(;F  lo  IIF(.  TM  470  IlFl. '  ) 
SFTmkn 
FNI1 

//DATA.  INPUT  nil  « _ _ 
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Program  IV 

Directive  Gain  Pattern  Calculation  and  Plotting,  3.17MHz 


CHhKt'  111  =  40. 
CHEkP I ? ) =?  7o . 
CHPkP (11=90. 


CHEkP (41=110, 
CHEpHlb  1  =  2*10. 
nn  .400  j=i  I .  i9 

Jl  NC=  I 

no  1 001  I = 1 , 1 H 1 

xval(  n=o. 


PHI «CMEKP ( J I 
PMIl«PHI-90. 

no  500  f= I .91 


*THS(  I  1  =  0. 

AOfll t 1=0.  , 

0  CONTINUE  '  , 


nil  2000  .12=  I  .  2 

THETAsO. 

nil  100  (  =  1.91 


TMETA=CHECK( 11 

FiNn  total  FiFLn,  fl.Aa  Is  factor  to  NORMai.I/F  to  (iAIN  nvFR  ISOTROPIC. 
_ nw.AF  (  theta  .PHI  I*FLF  (  THFTA  .PHl  1-11.62 


ATHAI  I  I =rHF  Ta 

AnB(n=(in  , 

TMFTA»THFTA*1  . 


CONT I NME 

wx  I  TF(  6,200  I  1  ATHA(  I  1  ,PHI  .Alim  I  1  .  (  >1 ,91  )  ' 

_ FORMAT!'  ' .51 F6.2. lX.Fh .2 .IX .FB.1.2X ) ) _ 

MRI TF ( 6.100  I  ~  .  » 

fohmati  '  '  I  . 

lP=91 

NrNjiANni n  ■■  ■■  ■  . . . - 

(X)  1001  1  =  1,91 


_  xVAi. ( ,iP )  «A  (iiAj  1  M.i  Inc 
.  TV  AUl  •il'  I  »  AIHH'  I  I  ■ 


■.iP  =  .iP  + Jl  M<; 

'linn  'III! 

'  IFT'ilNc.  "IhO  .  ly.ilMcs-  I 

Pill  .Pill  ,  I  Mo 

IMI'OI  .  K.ji.  l_rii|  =IMI|  -  if.ii. 
Cl'iN  f  I  MIIH . 

Ml,  I  I  I  I  10,1,.  I  I  jtWAI  III.  ,«A|  III, 
I  III, HA  I  I  '  I  I  M  I  I  /  . ,  •  ,  •  .  I  ,  .  I 
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■nn 3  tgit  1 I  I  c^LL  MV  tplti  iH.b.-oo..  jg.  .-ij.  ;3(i.  .  mn — 

IF  ( j  ,GT.  Ill  f;n  ro  7ino 

CALL  AHTPLTI  l«.F,-'-n..9o,.-T0.,30.,|Kl) 

Si  , 

2  1  00 

r  =  <»  ,bH 

03. 

X  =4 ,00 

ss , 

r.Ai.L  lfttfw  ( X  n  . ,  1  s ,  *  PMi  =  *,0.0. SI 

IS7. 

call  hiimrfR  (  a  .  V  .  ,  1  s  ,ph  iL  ,n  .0 , 1  ) 

SP  , 

X.Xa.95 

09. 

V=4.5H 

70. 

CALL  LFTT  FH  1  *  .  T  ,  .  1  F  .  MlFi;.  W'.O.O.AI 

71  . 

call  MFWPFMI X 1 

72. 

/•on 

CONTI  HI  IF 

73. 

call  finish 

7a. 

loop 

CONTINUE 

7S. 

stop 

7iS, 

"Twn - 

rr. — 

function  aF I THFTa.PhI 1 

7B, 

THFTAR  =  (THFTA/iHn.|AX.iFiS<J 

79. 

“PmIr* ( PO 1 / 1  ho . 1 »3 . 1 A  1 39 

c 

X 

.  r6  TAaF.FPPSAA  .  , 

HI. 

86TA4»*.»RFT A 

«7. 

rF  TaH.H .arfTa 

H3. 

SThETa^sIni THFTaR  i 

B4. 

SPHIiSiNIPMIh 1 

HS. 

AI=T«(SlN(RFT*4»STH«=T**SPHr  l/'SINtfleTA*STHFTA*SPHI  I  )  •  (  SI  N|  BE  T4ft*S  THP 

ITAACPHI  )/«;|N(nPTA»STMFTA»r.PH!  I  > _ _ _ _ _ 

"Trl  AHS(  IHF  FA-O.  I  .L?.  .X  .im~*HSMHFTA-lHO,  I  .LF.  .2  .UK.  AK^IHH 

ft-o.i  .lf..  ,2  .nn,  Afls(PMf-i«o. »  .lF.  .2)  i;n  in  loo 

r.O  TO  FOO 


HA. 

HI. 

HB  . 


S'). 

vn. 

91  . 


100 AFTl «A. ^ ^ ' ! 

GO  TO  300 

JOO  *FTlaABS(SlN(wFTA4«STMFTA«SPHn/SlN(qFTA«STHETA*SFHt  I  I 


IF  (  AnS(  TMpTA-o.  )  .lF.  .2  .I'R  .  ars  I  THF  T  A- I  H  0  ,  )  .lF.  .2  .OR.  AHSIPHI 
l-RO.I  .UF.  .?  .UK.  AHSI PMI -270, I  .lF.  .21  GO  TO  AOO 
r.O  TO  inn 


9?. 

93. 

9A. 


300 


AOO  aF T2sM . 

GO  TO  GOO 

_  Soo  aFt2 jahs isIn(rptah«sthFta»cpm| >/3|ni bet aaSTHFT a«cPH1 I  ) 
”^000  oraF T 1  sFo . » . A 34?‘*A»ALn(;( af r  i  i 
nRAFT?*20.*.A3A?9A«AUIl(;(  AFT2  I 
G _ 9R|TF|6.TnOI  THFTA.PHI  .AFT  I  .AFT2 


9S. 

9A. 

97. 


9H. 
99. 
1  DO  . 


ini . 
107. 
10>. 


Too 


FORMaTI*  •  .  '  THFTaa  *  ,FA  .2  . 'IlFG  PHI  =  '  ,Fb  ,2,  '  nFG 
12. <  S  FlFmFNTa'.FT  2.'  *1 

-A^.arsi^FIJ _ _ _ _ 

IF  (AOSIAFf)  .1,1.  I.F-IOI  i;(l  Ml  in  _ 


4  FLFHFNTa ' .F7. 


AF»20,*.A3429A*At,ll(,(AFI  » 

«nnAFTj  ♦hmap t p 

“*  ■■  ■  ■  "  . . 

10 

KHJIIMN _ _ 

FNn-'"  ■■  . 

HtNCnUN  FLF(  TMFTA  .f»M!  I 

_ «  M  I  .  A  I  Ml  f/\f  1  /  I  ,  M  (  /  ,  W,  )  p  I  F'i  f  ,  1  ) 

n/if  A  aPm  !/n,  .'i’o . ,  i  ui .  •  i  un  •  i  o  ,  *  ;>  7n, ,  A.m ,  ,  .  /  ,  a  im  \  f\  7n ,  ,  p  ,  ^ 

I  •  •  I  o ,  I  I  :*  p  *s ,  I  s ,  ,  I  /  p »)  ,,•!»  p  p  p ^  f  ^  ^  , 

1  AS  7  ,r,n  *-T ,  i.s  ,s,/.s  ,s  , 


.  I 


!  ,«  P.'i  .  HS,  ,  H  /  ,  S  , 

(I  A  I  A  I  I  /  f*  1  ,  ‘  » ,  /.  I  •  • , 

I  1’>  ,  1  ,  M»  ,  <  *  h  ,  <  A  ,  4  ,  /•  /  .  4  , 


W  ,  ‘i  ,  *iM,  ,  f. 


1  »>/.. 
lOS. 

J  . 

MU, 
MIH. 
1  •»>(  . 

I  Ml, 

in. 

I I 

1 1  ■». 

11/.. 

1 1  ** . 

I  M.  . 

in, 

1  I  M  . 

I  I  7  . 


I  ,  w , 


It -ft.  -4,. 10, -8. AS  ,-1^,00. -7. 19.-10.47,-11  .AH, -1.22,-4.88, -9, 78,  -  |  174. 
n. 99.-8.00,-!  1.97, -17. 81,-8.17,-8.88. -9.98,-14.54,-8. HA. -I  7. 98, -I  7  178. 
I.Ta.-9.09,-8 .57. -10. 5 5. -I  A. 88000. -9.  70,-1  .  . 90 , -I A . 8  1  ,  - 1 0 . 09 , -  7 . 44 .  1  7h , 


1-11.05,-14.87.-10.89,-14.87,-15.85,-11.15.-8 .47.-1 1.51,-14.77.-11,  177. 
154,-15.77,-18. 87, -12. 78. -9 .51, -12. 07. -14. 70. -12. 60, -16. 79. -17. 91,-  178. 
117. 57, -10. 77. -17. HO. -14, 84, -17, 90, -18. 01, -19. 11. -14. 99, -12, 21/  17q. 


nn  700.  I  =1 ,74  , 

ml  400  j=  1 . 7 


i«ti( j.i  1 
CONTINIIF 
CONT  IMIIF 


IMI  500  1.75.77 
no  AOO  Jxl ,7 
I T.t-74 


ep ( J . I |.T2( J.ITI 
CONTlNIlF 
CONTI MUF 


on  loo  I. 

17. 1 

IF  ( A8HI  ( 1 1  .r.p .  PH  1 1  r.n  to  looo 


I  I  .  I  I  8  I 
I)  CONTINIIF 

1  no  700  J.2.78 _ 


J7.J  ' 

ifmthetaui  .gt,  rHFTA)  r.n  rn  iioo 

ji.ji+i 


WR I TF( ft, 7000  1  ( ( FF( I .J ) , 1  =  1 ,7 ( . J»l , 77  I 

7000  FORMAT^  ',771'  •  ,  7  (  F  7  .7 , 7X  I  ,  /  I  I 
Jl  =  l  _ 


PH  11  r.n  TO  1000 


SET  J7  FUR  lInfaK  F  X  TR  amiil  a  1 1  Iim  mil  Til  *rO' IIFURFFS. 

.I7.7ft 

J1.J7-I  ^ — 

7  FR  act  .(THFTA-aTHFTAUI  (  l/'(AlHFrA(J7l-ArMF|A|JI|li  ' 

_ FRAC8  =  |HH|-8PmI  111!  1/  MPlil  I  I  7  1-APHl  n  1)  1 

IF  07  .pij.  hI  1  7^  I  ■  "  . 

EF1*FHaCH*(FF(I7.ji|-FI-(|i.,ii||*(.|.|Ii,j|) 

FF7.8RArP«(  l^|_U  ,  17  l-ll-  (  I  I  .,1  7  I  I  HP  (  I  1  .  Ml 
FLF«FMACT»7FF7-Ff- I'liFFI  ■  " 

WR  I  TF  1  8,70ll  I  FI.F  .A  Mir  1  A  I  .M  I  ,A  III!  |  A  |  II  |  ,  A|>|l  I  |  |  I  1.«l.|l|  (  |  |  | 
IUHMA|(_I  O^I^II-MIMIAI  I  A<  |IIM^'  ,1  8.7.  '  Ull  •  .  4  I  I  X  ,  *  A.  7  |  ) 

8F  niMN 
IMIl 

•.I  IllMl  Mil  I  PI  I  Afl  I  I'  I  I  I  l-M  IH  •  ^l|  Mill  ,  M  I  M  ,  IM  A  <  •  yM  I  1 1  ,  y  M/i  »  ,  |  •,  ) 
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CIIMMIIM  /'BLnCxi  /xval  (  1  «  1  I  .VV'AL  M  h  1  I  •  I  Vl^dliJ  )  .  1  XIICIMH  ) 
niMFNslilM  X  vvw  (  .  yuu\j  ( 

INTEliFK**  VSCAI.F ,  xsr.Ai  t-  _  _  _ _ _ 

“  “  NUMsi'fiF 'll  I V  fs  ImaT";  mm  x  axis  "  '  • .  . 

C  MtIMn  s  »  IlF  nlV/ISlMMS  MM  Y  AXIS 

r,  XM|N  A  VALIiF  MF  X  A[  MHKIIM _ _ 

r.  XMAX  »  VALMF  MT  X  at  FMM  MF  AXIS 
f.  YM|N  =  vai.mF  of  Y  at  miiIIUm 

0 _ YMaX  =  VaLMF  MF  V  at  TMP  IIF_Y_AXI_S _ _ _ 

r,  nfnts  =  n  riF  I'mImts  Tm'hf  plmifm.  must  hf  i.Fss  imam  «T 

CALL  FLTTYF(SA62.h .T7 I 

_ CALL  STaKT _ J _ 

CALL  HlOTIO. 0,0.0, XI 
CALL  FlOT 10.0. 1 .0,2  I 

_ CALL  FLOT10.0.7-.Q  .XI _ 

CALL  FlUT lO.O.H .0.2  I 
CALL  FlUTI  I  I  .0 ,0 .0  .X  I 

CALL  FlUT  I  10.0,0.0.2) _ ^ _ 

CALL  FlOTI I .0,0.0. XI 
CAI^L  FlOTI  0.0, 0.0. 2) 

_ CALL  MFWFFNI  I  I _ ^ _ ; _ »■' 

T,  nfcFINF  mFW  mmIijIm  fiif  plot  axis 
CALL  FlOTIP.22  ,  I  .HH  ,-XI 

_ CALL  FacTIIWI  I  .01 _ _ ' 

r,  iiRaxJ  AXIS 

ENTRY  NXTPLT(MIIM,NIIMn,XM|N,XMAX,YMIN,YMAX,NFNTS) 

CALL  HFCTI  0 .0,0  .0 .7.0  ,S  .o  I 

7  flKAW  Tic  maHks  mm  axis.  CSI7F:I>I  visiliN  SI7F  IM  IMCHFS  X  AXIS. 

r.  osiZ6»ni VISION  sUf  1m  imchfs  y  axis 

_ CSIZFaT.O/PLOATINI  IM  I _  .  . _ 

nSl2F*&.02FLOAT(NIIM0l  '  ,  , 

XI >0.0  '  •  ,  ‘ 

_ Y  1>0.0 _ _ _  ' _  •  _ 

T?>Y1*.0S 

XPsXlA.OS 

_ nil  Aoo  KA 1 ,2  _ 

no  200  j>i  ,NiiM 
•XHASF  =  FliIaT  (  NIIM-\j  I  ACS!  ZE 

_ X  aXRASF _ ■ 

CALL  FlDT I  X .Y 1 . XI 
CALL  FlOTIX.YP.PI 

200  cuntinmf _ 

CALL  FlOTIXI  .Y1  ,XI 

CALL  FlOT|X2,Y| ,21  i 


YMASFsFLIIAI  (MMMM-.I  I  *MS|  ft- 
CALL  Flu  1 1 X 1  .YHASF ,X 1 
call  F|.|IT|X2.YHASF.2I 

» . 

’  a 

• 

^oi) 

cuntinmf 

X) >7.0 

Yl«5.n 

'  yp.Yi-.os 

400 

rflNTlMiiH 

‘ 

r  FUT  SCALP  I  IN  Y  AXIS 


0.  UNITS  B  UNITS  Fill  lllv  MN  Y  AXIS 

_ IIN1  Ts«(  (  YM  aX-YM  In  I  /I  I  MA  I  (  MMHIl  1  I 

YSCAI  F  >  A  I  N  I  I  AMSI  Ym  I  M  I  I  .S  I 
II  lYMlM  .11.  I  Y’.f  Al  I  Y'.,  A,  I  ”1"  I  I 
_ _ !.li;.=  !lMMLM.l _ _ 


TwTT 

1»X. 

I  H4  . 

jws. 

I  FA. 
1F7  ■ 

I  FH. 
|H4. 
JVM. 
IVI. 
IV?. 

I  vx. 

1  V4. 
IVS. 
IVA. 
IV7. 
IVfl. 
|VV. 

?i>n. 
?nl  . 

?nx. 

?n*. 

?I>S. 

?06, 

7>n7. 

?nw. 

20V. 

?in. 

?l  1. 
?!?. 
?1S. 


?1 « . 
?IV. 

??n. 
?7I  . 
???. 


2?x. 


2 


.  I  f 


777. 

77H. 

7?V. 

?xn. 


7X1  . 
7*?. 
7XX. 


?XA. 


2  IS., 
?XA. 

2  •‘PI  . 


2  <'l  . 

>  I 

•  • 


lh(y<;CAI.f:  .I.P.  m  <  =  -.^07 

CAi.L  I  miimhk  (  X  ,1  iLiir  . .  1  s .  Ysr  Ai>  ,1) .  n  I 

Y<;CAi_F  =  Y'^r4l.(^*AlMTIiiM|TS*.‘>  • 

400  cilNT  I  MMF  ■■■ 

XT=-.5 

_ YT»1  .b _ _ _ _ 

CAI.L  LFTTPK(XT,YT..l‘i.l  YIICI'  .•>').  O  1  A  I 

r  PUT  scalf<^  um  x  axi*; 

C  IINITX  UNIT";  PFK  1)1  W  DM  X  AXIS 


LARfc  L  X  AXIS  SCALPS 
I  T  1  =NIIM+  1 

no  Ton  I  =  1 . 1 T 1 . s 


IIL0C  =  CS|  7p*(  I  -  I  I-  .  l?S 

IF!  xscalp  .lt.  0  I  iiLOCAiiLnc-.  IS 

CaCl  iMIlMwa  IIILOC,  -  .?A  .  .  1  S.XSCALF  .0 ,0  I 


XSCALF»*SCALP*aINI(IIW|  IX*.S  I  ax 
cnNT  1  Ml  IP 
XT  =  I  .0 


call  LFTTFHI X  T.YT. . is, ?XHCn,0.0,X2) 

IIMiTS  PPW  niv  OM  Y  AXIS _ 


Plot  oata 

ENTBY  PlIMFIMPNTSI 
no  Ano  I  «  t .npmts 


IFIXVAli  I  I  . 
IF(XVal(I|  . 
IF( YVALI  I  1  . 


IF  (YVALI I  I  .  ,  ^  . 

CONTINIIF 

WR  I  TF(  A  .  I  on  I  I  (  XVAI.  (  I  I  .YVALI  I  I  .  I  =1  .  I  HI  I 


YMAX I  YVALI I IaYMAX 


CALL  PLOTIX.Y.SI 

NPNTS  *  A  OF  POINTS  TO  RF  PLOTFO 
no  Son  IsP.NPMTS _ 


XsIXVALI  I  I-XM|NI/(IIN[TX/CSI/FI 
Y»  I  YVal  I  I  1-YmIn  )  /  (HMl  TS/l)SI  7P  I 
CALL  Pl.UTlX  ,Y  ,2) 
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ENn 

HLOCK 


n«TA 


CIIMMIIN  /Ml.dCK  1  /X  V*l.  (  1  M  t  I  ,  VWAl".  )  I  hTi  ,  I  YhCHIA  I  ,  I  XHCIKH  I 
DATA  I  YHcn/ ‘HFaT  '  . 'E«  t;i,*AIM  •.MIIHI'/ 

MATA  I  xncn/ '  '  « ' _ *  I ‘The  I*.' A  ' .  *  idec,  ' . 'he^s'  , 


■XlX. 
^IIA  , 
X(TS. 


I 


1/ 

ENn 

//nATA.FT^7Fooi  PIT  »dL*REPJ  MFN  .llfl4A40  .KJC.LIH.DISPstNFW.ItFFFt  . 
//  HSN=MFN.(IH/./.'Tn.KJC.F(.nT  .riflTPHT  .SPACEMCYC.  (1.1)  ,«LSF  )  ,  ^ 

//  tlCHs  (RECFMsFR  .I.KFCLsHo  ,HLKS  I  /  F  =  1  IPO  ) 

//OATA.  INPUT  nil  • 


^(14. 
>(i7, 
XflH  . 

Vtl'i. 

^in. 

’ll. 


?i  7. 
’13. 
314. 
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//SCH  fcXHC  t'(.M=I  PFMft  li  ^  T" 

//sr.«  nn  ni  sf'=  ( nLi'>.nFLFTP  i .  vtn.=«FF=u(''.o'>i ,  • 

//  l)flN  =  MHN  .KJC.  PI.  (IT  .mil  Ml  IT _ [ _ 

//  i-XPC  l-XCI*  "  ■ 

/•JH  SFRV/ 1  rF=t)^FFK 

/».)F  FIII.I.SK  I  Fs _ I _ 

✓/SVFIN  nn  * 

CdMMUN  /Bl.nttK  1 /X  VAI.  (  IB  I  1  .yVALI  l«I  I  .  !  VHCn  14  (  .  I  XHCnjH  I 
r.  4  CHABAcTFHS  FFR  IMTFr.FW  in  I  mTFr.FW  AWRAYB,  SFF  NtflCK  DATA  SDMpRfir,  I 
nl  MFNBlIIH  ATHA(  S  I  I  ,  Aumm  tVCHFCK(  S  l  .CHFKP  (  S  I 
CHFck (I  I *n . 

_ CHFck  (  ?  1=  _ [ _ 

CHECKMIiFF. 

CHECK(4  )»S7,S 

_ CH6ck(5  1 1=90  . _ ] _ ^ _ 

CHeKR|M=40. 
dHEKF { ? 1 =?  7o . 

_ chpkp  ( 1  i-ggn. 

CHEkH ( 4  I = 1 ^0 . 

CHfKt»(l)  Ia^SO.  ■  • 

_ nd  4oft  .i»i  _ ] _ I _ _ ■ 

.|IMC=I 

nn  jon'x  1*1 ,  iM| 

_ XVAi.  I  I  l»n. 

100^  YVAL(  I  )a0.  ]  ■  !  — — — — — 

; _ PH1»El0aT(  J1  l«l<7.»<70» _ _ ^ _ 

PHlL«PH|-qo. 

_ nn  500  la  1 .9t _ _ ; _ _ 

atha(  I  )■<»,  ^  ^  ' 

Ann(  n«o.  '  • 

5dO  CONTINUE  ■ _ ■  _ ^ _ _ . _ _ 

n»l  ?000  J?x|<7 
THETA. O. 

_ rxi  ion  1.1, qi 

"T  thSta.chfck  ( 1 1 

C  Finn  TiiTBU  FiFLn.  q.H?  is  facthb  Tn  nobmalITF  Tn  gain  over  isnTRnpic 
_ .  oh.af  I  thfta ,fhi  i»pi.fi  tmft a,fhi  i-q.wg _ 

ATHAI  I  I .THFTA 

AnR( 1 1 .nn 

_ THFTA. THFTa*I  , _ 

100  CUNT  I N1IF  : 

WHlTFI^tPOOl  (aTHaI  I  I  ,>M|  .Ann!  I  ».  t»l  ,911 
700_  FIIHMATC  •  ,‘>ltF.?,IX«FF.7,l*.F«.R.7xil. 

mrTtfio.^ooi  ■  ■  •' 

qoo  FIIHMAT  (  '  '  I 

_ JF.9  1  _ _ _ 

NENn.NHT-fl  ^  ^ 

nn  tool  i.i,9i 

_ XVALI  JF  I  »ATll»j.!  lAjlNC 

yvalijf I » Ann rn  '  .  .  .  _ 

jF..|Pa  Jl  NC 

1001  ClINTiMIIF  _ 

I  F 1 .11  Nr  ,FiJ,  II  .il  Mr  .-I 
Fill  .FMl  A|  Hn 

_ IF  I  Phi  .i.t.  ^fo.i  mmi  =i'mi- 

Xorin  rllMlfMin- 

Mil  I  I  I-  I  >1  ,  I  IIIIY.  I  I  XUAI.  I  I  I  .  VUAI  I  I  I  ,  I  :  I  .  I  II  n 
lOOA _ FIIHMAJII  I  (  III  F7  •  >  .F  7  .X  I  . /•  'll 
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1BPTAH««;thp 
•  flR,  flRSIPH 


call  lfttfri  x.v.  .1  r  . 'nFr..  ut.o.o.M 

CALL  NFWPFM(RI 
I  CUNT  I  All  IF _ 


CALL  finish  - - - - - - - - 

P  CONTINIIF  ,  , 

srnt* 


Nn  : - - - - — 

functiun  AF ( thfta .phi  ) 

THFTaKs ( THPTA/t Mo. ) wx. 141  so 


PHIRalPHI/lMo.Ias.lAlSO 
RFTA=2*PI*(F/CI*HS/2.  AS=n|sTANCF  HFlwFFN  RAniATHRS. 

B6TA»4.SS96ni 

""i  BFTAAsA.ARFTA  - - — 

BFTARaH , *BFTA 
STHFTAiSINI THFTAH 1 


SPHI aSiNIPHIM ) 

CHH|.cOS(PH(H I  * 

>/RINlRFTA«STHFtA»SPHt  I  ) »  (  S I  N  I  BF  T  AH  «  S  THC 

ITaacPHI  i/sIn(bfta*sthfta»cphi  )  I  — *- 

IFI  ABS(  THFTa-o.  1  .lF.  .?  .flP.  ARSITMFTa-180.  I  .LF.  .2  .fIR.  ARSfPH 
U-0.)  .LF.  .?  .OR.  .AHS(PH|-IRO.  )  .lF,  .2)  F.fl  TO  100 

r.O  TO  200  ^  ^ - 

aFT«=a.  •  , 

on  Til  son  ■  ‘  ' 


AFT  1 s ABS( SiN ( RF  T A4*STHFTA»SPHI  I /S | N ( RF T AASTHF T A»SPH |  1  ) 

1  IF(ars(TMFTA-0.)  .lF.  .2  .flH.  arsITHFTA-IHO.  )  .lF.  .2  .PR,  AHSlPHl 

- L-lRfl  .I.F.  .?  .PH.  AHS(PH|-2T0.  I  .LF,  .2)  i;i  I  TP  400 

GO  Tl)  boo  :  - — - — — 

AFTPaR.  ' 

_ liO-  TP  AO  o _ _ 

AF  r2«4HS(SlN(RFrAM*SrHFTA*CPHI  I/'SINIRFTA»STHFTA*CPHI  )  I  - - 

P‘'F  .4A4204  TP  CPNVFMT  NATiiHAL  LPi;  IP  HASF  10 
PR AF T I a2n . «.4 T4P04AALPG ( AFT  I  I _ 

nBAFTP  =  P0.<>.4X4?94*ALPG(  AF  T2  I  ^ ^ - - 

PRITF(6.700I  thfta. PHI .PHAFTI .PRAFT? 

•~"r;  '■  ^^'**'^* T I  *  '  t '  THF.I  «;JaE.Rj?t.*.DEG. _ PHlJLLi£?2.Lej_il'Fr.  _ ft_Fl  Fhfnt.  ■ ,  f 

I?,'  PH  H  FlFMFMIs' .F7.?.'  PH'  I  *  - ^r.i  r  i.  ,,r 

aft>ahsi  aft  I  . 

_ If _ (ars(aft)  .i.t.  i.f-101  (;p  tp  in 


_ Lf _ (ars(aft)  .i.t.  I  .F-ioi  f;p  rn  m 

AF«P0.«.4S4294ttALP(;(AF  r  I  ~  - - 

aF  xPRaF  TI  aPHAF  T  7 
RETURN 


10  .  aFs-WO.  ■■  ■■■■ - -  — . 

RFTHRN 

_ ^fnp 


function  FlFITHpta.PHI I 

PI  MFNsIPN  FF(  'I,  17  I  .AI'HI  no  I  .AlllFI  A  (  1  7  I  ,  I  I  I  7  ,  1  ,  I  ,  I  ^  I  ,  ,  ,  ,  , 

!!,*1A  AP"|/|,,  ,4„.,‘(0..|  . . 711.,  1^(1.,  wo./,  A  I  HI  lA/n 

I  ..7.S.S..  7.S.  I„..,;,.S.  7.  S., ,  17.S.  IL.  .  w.s 

I  .411  4  ,4  7  .Ml  I,  ^,1,  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

- f-9^  <  «  S_«  M  S  ,  ,  H  /  ^  S  a  «IM  -  y 
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nAi'A  T 1/ '*1  .K7 .  ^ .  /  A.'^ .  .'A.  /i.T f. ^h, 

I  A4  ,  A.HA,  A  .'5?.  .  A  ,  A.i.1)  ,  A.SO,  ,  A  .£)?,  A.<<  7  011(1.  A  ,  ?H  ,  A  ,  . 

lA.?4.A.A4.A.7<..A.HH,A.S7«A.|q.?.44.A.O?.A.|S.A.h().A.H?.j.VA.?.‘<A,  ^ 
1  .  AS  ,  ?T7A(in,  ?.S  A,  A.S4  ,  A  .7^>.  A  .»H  ,  ?..\H  ,  AV  .  ?  .(,((  .  A./,(l,  A.  f,;  A.  74  .  ? 

l.AA,^.q<,.^.l?,^.A4.A,2A.A./,^,A.^.H.1.44.1.S?.l,74,^,07,A.()A,A.?4,A. 

lSAn,l.S7.  l.ns.  I.Al  ,  I  .70.P.H1  .A.OI,  A.AA.I  .«q..AP,.l«2.1.?4.?.SA,7.7S 

I  .  A.  IS  .  .AS.  -  .OA  .  .?H  .  .BA'.P  ;?H  .?  .44  .2.9  I  ,  -  .OA  .  -  .  7  1  .  -  .A2  .  .33  ,  1  ,qHV?.  I  a 

1.2.64.-.74.-1.4I,-.94,-.?4.1.AA,I.7A.2.34.-1.S1,-2.19,-I.7A,-.H7,I 
1.A.I.AA,2..-?.A7,-A.OA.-2.AS.-I.SA.0.9A.0.9A.I.A2.-3.34.-3,9A,-3.4 
lA,-P.AA,.SS..4A.1.21.-4.42.-4.9A.-4.47.-A.lH,n.lS.-.n4,.77,-S.44.- 
15  .9  7 .-5.A9 .-4. 1 A,-.27 .-.S7 . .29.-7.O.-7.0A ,-A.«2,-S.  1 7,-.7n, -I  .  1400 
1.-.22.-H.SA.-H.I9.-B.I9.-4.A4.-I.14.-1.7A. - . 7  A . - I  0 . 22 . -9 . 30,-9.700 
l.-7.AA,-l.S9,-2.AS.-l.A4,-tP.0S,-l0.A(.,-11.34,-9.07,-2.06,-2.29,-l 
1.94,-13.90,-1 |.A,-I A.t0.-10.Aft,-2.SA,-3.ftS,-2.S7,-IS.4A,-12.OH,-I4 
I  .(1700.  -1  2 .47  ,- 3. 1  ,-4 .34,-3  .3?  .-lA.27.-12.A5.-IA.4B.-lt.40,-3.-AB.-S 
I.04,-3.H9,-14.|9,-|A.0S,-17-A2,-IA.39,-4.32,-S.7A,-4,S9,-1S.S4,-1A 
1.32,-l(l.ii,-iH.lA,-S.04,-A.S,-S.32,-l4.7(l,-13.S4.-lfl.nA,-19,37,-S. 
1(1A0000,-7 .2A.-A.n4.-14.1S.-t3.79,-l7.«.-19.AA,-A.81,-fl.07,-A.9I,-l 
13, (H, -14. IS, -I/. Al. -19, S, -7. 93, -R. 94. -7. HI, -13. 7A. -14,71,-17. 7, -19 
1.330,-9.30.-9.9S,-M.fl7/ 

n*T>l  T2/-14. 23, -1  S. AS, -IR. 25, -19. S«.-l  1.05,-11  .22,-10. 21, -IS. 57,-1 
1  7,  A4',- 19. 69,-20.7  2.-1  3 ,52.-1  3.  1  3 .  -  1  2 . 22  , -1  9  .  1 2  , -21 . 1 9  , -2  3 . 37  , -?4  .  1 
17,-1H.01 ,-l 7. 03, -16. 23. 7 *0.0/ 

OATA  T3/3.H3,3.H3.A.(10,3.flS,3.7S.3.H4,3,6(1.3.«P,3.5(l,3.76,3.47,3.6 
l(»,'3.32000,3.s7;'^liA.3.44,2  i47, 3. 27.2,76,3. 0  7, 2. 52, 2.  (13. 2. 25,?.  5600 
1 , 1 .96,2.24, 1 .64, I .H9 ,1  .29 ,1  .49; .92, 1 . OS , 0. 52 , 0 . 56, . I  0 . .03 , -. 36 , 5 
15.-.H4.-|.i(1.-i.3S.-|  .HS,- I  ,(19.-2.56,-2.45.-3.30.-3.05.-4.06,-3.66 
1  ,-4.84.-4.36,-5.61  , -S. OH. -6. 37. -S. (15,-7.09, -6. >0,-7. 76, -7. 64,-8.40 
1 .-H.69,-9.01 ,-9,89.-9.66.-11 .33,-10.43.-13. 13.-11.49,-15.66,-13.74 
I .-20. 2200, -I  7. 07,0, .0./ 

J__  ;  ;  ^  , 

I(A(  tAMI  .lT.  10. I  PHI»PM|4360.  <■  , 

no  300  1  »1  .33 _ /  _ '  _ 

no  400  J«l,4  I 

EP( j.i  )>Tiij.n 
400  CONTImiif 

Ef!(4,  n.T3(  1 .11  ~  ^ 

no  450  J»5,7  .  ' 

EF(J*l.ll»Ti(j.ll _ 

4^0  CONT INIIF 

EF(9, I  1.73(2,1 1 
300  CONTI  Ml  IF 

no  500  1.34.37  :  [  -  ^ 

no  600  Ja  I  ,4 
_  IT»I-33 

FF(J,I  |=T2(J,IT1  '  »  . 

600  r.ONTIMIIF  v»  , 

_ FF (5,  I  ).T3(  I . I  I _ 

no  650  .1.5,7 
(:K.I+I  .  I  I  .7  2(  .1,  I  I  ) 

650  rllNT  INIIF 

l->  (9,  i  I  .7"3(  pTTT  "  "  • 

500  rONTlMIIF 

r  WU  I  (F  I  A.piliidl  I  I  I  F  I  (  ,  I  I  .  I  •  I  .<1  I  .  .1  I  .  1  7  I 

r.  2000  i-OHMAri'  'V^Ti'  '  .vift.p  .33 1 ,/ 1 1 
.11  •! 

no  Inn  I  -  ,  I II 

M-l  . 

I  I  I  Al'll  I  III  ,1,1,  I’M  II  l.ll  III  I  IMMi 

I  I  -  I  I )  I 


rmi  riiNn  wTF 
lono  nf)  ?on 

_ - _ 

IP(  AlHf-'TA  I  J  I  .(;i.  '|KHT4l~i;n  III  lino 

J 1  =  J I  ♦  I 

JOO  CllNTlMMF 


Jl = J^-1 

1  I  on F  M  ACT  =  (  IH  PI  A-ATH  F  1*  1^1  j  j  .1  /  1  *  1 IJ  PTA  I  J 1 1  -  A  i  ^  J  !  L  • _ 

IF  111  .FO.  1^1  WK  1  T F  I  A . Hoo  i  PH  I 

F»ACF= (  FHI -aPH  1  (  1  I  I  (  APHM  n  I-APHI  I  I  I)  ) 

_ IF  I  n  .Fu.  10 1  n=i _ 

EFl=FHACF*<FF(I1,Jil-FFIM.Jin*FFIIl.Jll 
.  EF?=FRacP*I FF (  I ^ . Jl  l-FF ( I  I • I  l+FF I  I  1 , Jl  I 

_ FLFgFRACTAIFFg-FFi  I^FFI _ 

C  WK I TF I h, 7on I  FLF.aTHFIai j^l .aTHFIaiji ).aph| ( m , APMl (Ml 

r.  700  FUKMaTI'  '.‘FlFMFMTAl  FACinR  =  *,FA.2.i  DB’ I  ' 
Ron  FORMAT!  '  I, Fin. 1,1  ph[  miit  IM  Ramof  i‘o  liFf.  TO  370  nEU') 
RFTMRM 
ENH 

_ nllnllQIlTI  MF  AMTPLT  (  MUM  .MIIMO  UMAX  ,YM|N  ,YMAX  .MPNT^I _ 

cOmUhn  /Bi.ncK  1  /  X  vAi.  (  1 H 1  I  .yvAi.  ( I  HI  i .  /hco ( *  ) ,  i  xhcihh  ) 
niMENFllIM  Xvww  (  3A  I  ,YVVy  (  3A  ) 

IntfCiFRaa  Y'lrAi  f ,x«;cai.f 

T;  ^jiimak  o»^  III VI shins  om  x  axis  ‘  .  ; 

C'  NUMO  «  *  OF  nlVlSlONS  ON  V  AXIS  .  * 

r.  XMIN  .  ValOF  of  X  AT  ORir.lN _ _ * 

7.  XMaX  «  VaI.OF  DT  X  AT  FMO  OF  AXIS 
f,  VMIN  »  'IalIIF  of  Y  at  no  If,  I  at 

r,  ymax  «  VAi.np  OF  y  at  top  of  y  axis _ 

"c  NFNTS  »  «  OF  POINTS  TO  SF  RUOTFO.  MOST  BE  LESS  THAj(  8> 

CAI.U  PlTTYI»(*6A2,A,37I 

_ CALL  start  _ 

call  Plot  10.0,0.0,3 1 
CALL  PLOTI 0.0, 1 ,0 .?) 

_ CALI-  Plot  (0.0,7.0,31 _ [ _ [ _ 

CALL  PLnT(n.O,H.o,?l 
CALL  PLOT! I  I .0,0.0.31 

_ £L*LJ _ PLOT!  in.O.P.o.gl _ _ 

CALL  Plot (■  i  .0 .0 .n,  ii 
CALL  Plot  I o  .0,0.0 .? I 

_ call  mfrpfm! 1  ) _ 

C  OEFINF  NFW  nKIf.lN.FOR  PLOT  AXIS 
CALL  Plot  12.22. 1 .RH. -3 1 

- C  ALL  F  AC  inn  11. .01 _ ^ _ 

r.  ORAM  AXIS  „ 

entry  NXTPLI!  mom  .NIIMO  .XMIM.XMAX.YMIN.^IMAX^MI'MISI 

_ CALL  RFCT ( 0 .0 , 0 . o , 7 . n. s .0  1 _ 

r,  DRAW  rlr  marks  om  A*Is.  rsl;F-0|  VIMOM  SI7F  IM  INCHES  X  AXIS, 

f.  IIS  I  /E  kO  I  V I  S I  OM  sl/l-  IM  IMtIII',  y  AXIS 

_ £5JZf  “.'.tO'lFiojAT^lfJOM  I _ 

ns  i  ^F=‘>  .n/Fi-iiA  r  I  NOMo  I 
XI. 0.0 

_ T,Li!5f.n _ _ 

Y?.Y 1 A.os  ■  ■ 

XP.X I ♦ .OS 
MU  Alio  K  _ 

IHI  pfjn  .i~  I  ,  Ml  IM 
*  MA'.t  -  t  I  MA  M  MUM-  I  »  •«  M  /  I 
X  .  XttA'.l 


13 


-V 


CALL  1  X  .  Y  1  ,  ^  » 

CALL  Hi  (J IJ  X  ,  y  ?  ,  ?  ) 

7no 

CtlNT  I  wilf^ 

CALL  Hijl]  (  X  I  ,  Y  I  ,  ^  ) 

CALL  Ht.ll  r  (  X?  ,Y  1  ,?  1 

rxj  ^00  j=  1 . Ml iMi> 

yMASP=i-L‘'AT(MiiMn-.M««>«;f;F 
call  ►*Lf)T  (  X  1  ,  YhasP  .  ^  ) 

CALL  HLllT  1  X?  ,  YhaSF  •?  ) 

-ton 

i 

Cl  INTI  Ml  IF 

X  1  =  7  .n 

Y 1 =S  .  n 

1 

X?=Xl-.0S 

Y?=Yl-.OS 

1  400 

CUNT  INIIF 

c 

PUT  SCAI.F  IIM  Y  Axis 

r. 

UNITS  =  UNITS  PFM  nlv  ON  Y  AXIS 

UN|TS=((YMAX-YMlM|/FLnATI  Ml  IMO  1  1 

Y SCALF =A 1 MT 1  AHS (  VM 1 M I +  .S  ) 

IF(YM|N  .lT.  01  YsCALF  =  VSCAI.F«(0-I  1 

IT?=NIIMn+t 

rxi  Ono  1  =  1 .  1  T? 

(ILflCsnsI  ;F4(  I-l  l-.OA 

Xa-./.ft 

IFIYSCALF  .(iF.  01  X  =  -,X07 

call  INIIMBU  (  X  .IILOC.  .  IS  .YSCALF  ,0w0  1 

YSCALF»YSCAI.  F+aINT  (UNI  TSa.s'  I 

^loo 

CUNTINI'F 

XT=-,b 

YT«t .S 

•  CAI.L  LFTTF«(XT.YT..Ib.IYBCn.<»o.o,IM  ^ 

<■,  PUT  SCALPS  ON  X  AXIS 

r. 

PUT  SCALPS  (IN  X  AXIS 

c 

UNITX  =  UNITS  PPK  uli;  (IN  X  AXIS 

IINITXal  (  XMaX-XM|N)/FLUAT(MUM|  1 

XSCAl.Fi  Al  NT  (  ARSI  XM  (  N  1  +  .S  I 

777. 

IPIXMIM  .LT.  Ill  XSCALF  =  XSCALF«(0-1  1 

77b. 

c 

LAhFi.  X  AXIS  SCALPS 

77Q. 

Tin 

I  T1  =NIIM*  1 

7Hn. 

nu  7no  1  =  1 . 1  T  1  .  s 

?H  ]  , 

ULUcacsi ;f«( 1-1 IPS 

,  7«7. 

IP(  xscAi.F  .lT.  n  1  iii.uCaiH.nc-.  IS 

7HX. 

r,*uL  I miimhh  ( Munr.-  .  .1*1  .xscalf.o.o  ) 

««iCAI.F«XSCAI.F  +  AlNT(tlM|  TXA.S  I  *X 
Inn  riiNTlNifF _ _  _ 

xi=i.'(r  ■■  " ,  . . 

VT»-.AA 

CAI.I.  IF  TTFH  (  X  T  .  Y  I  .  XHriijii.i).  _ 


_* _ _ 

V  . 


IIWITS  «  liwfls  Mf'k  lllv  (IM  Y  AXIS 
FI.IIT  flATA 


F  M  T  H  Y_ .Fl.J  M FJ '‘FJ'' T  •;  L. 


i«i  Aon  i  =  1 
IF(XVA(.(  1  1 
IP  (  XVAI  (1  1 

.mPmIs 

.If.  X  M  1  M  1 
.1.1.  XMA  X  1 

X  VAI 

YVAI 

(  1  ) 

.1  M 

7  V  M  1  M 

2  XMA  X 

IP i  yvai  ( I ) 

,  r  1  ,  y  M  1  N  I 

YW  A| 

1  1  I 

-  V  M  1  M 

II  (YVai  1  1  1 

.1.  1  .  YMA»  1 

VW  Al 

(  1  1 

•  vriA  X 

400 

(  IIM  r  1  ^M0 

WM  1  1  1-  I  ^  t  1  OM  1  M  X  ^>AI  (  1  I 
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;  Y  =  (YVAL(  imI  TS/lls  I  ;  F  1 

CAI.L  ►'LUT  {  X  ,  Y  .  ^  ) 

r _ Nt^NTS  =  «  HP  HlllMTS  TH  HP  PlIlTfn _ _ _ 

Dll  5on 

x  =  (XVAL(ll-XMlMi/M'Mii«/r';l7Fi 

_ Y-IYVaiiI  I-YmImi/iumii-^/iisI  1  F  ) _ 

CALL  I'LUTIX.Y.?! 

Sfin  C.IINTiMIIF 

^‘LCIT  RICF  PaFFK  data  _ _ _ _ _ _ _ _ 

Data 

r  in.,-lh.?,-|F.'Y.-lX.».-in.X.-A.H.-A,lH,F.2(l.0..A.A.‘J.M,ll.^,14.X,IS 

f. _ l.Wt|7./.?n..?X.A.?S.H«gH.A.Yl.»,X4..<A.l.lW.n  ,  .4ii.3,4H.A.SX.F/ 

~y.  DATA  Y  VV/M/7 . 1  1  .-1  .X*!  ^  .HH,  I  «i.6  1  ,  IS  .M  ,  1  1  .ftl  ,  -  -  1  ,H4  ,  1  1  onn. 

r  I-?**  1  I  ,  14. >5  I  I  |  ,H,  1  I  ,-7.74 

r.  I  .  -7.H4 , 7  .  1  1  .A  .A1  .ftl  .-I  ?  .I<4  .3 .7  I  .Y..M  .-.HQ,  -  I  7  .A«r,  -1  .34 , 7  .  1  1  .  I  7. 1 

T.  I  I  .  I  7.  HI  ,  I  A.  1  1  ,  1  I  .A  I  .  1  .hi  / 

r  n(i  100?  I  =  1  ,  7A 

r.  X=(XVVy(ll-XMlM|/(IIMITX/C4I7FI 

“c  Y  =  (YVVV(  1  I-YM|M|/(IIMITS/I1S17F) 

r.  CAI.L  ■'iYMHnLi  * ,  V ,  .0*;  ,*,o-.o n 

r.  loop*  CIIMTlNUF _ _ _ ' _ _ 

~f.  whTFfi  o,  1  ooiTT7uvv7TTT7yvuw(  1  I ,  I  =  1 .  7h  i 

dFTIIMN 

FNn 

HI.IICK  DATA  . 

ClIMMIJN  /nunCK  1/XVAI.  (  IHI  I  .VVAL(  l«t  )  .  lYHCHI  A  I  .  F  XHCTIH  » 

__  _ .r!*TA_ I XHcn/ ' mpat  ■ ,  ■  fh  (;» ,  ■  a f m  j_,j  1 011)3.^; _ 

"ilA  (  A'lSdicii/'  "  iilF  r*  .  •  A  '  I*  ,  '  (  IlFi;  '  .  'HF-FS  •  .  •  I  ' 

I  / 

_ Fr.M _ 

//data  .FTiVfiOO  1  nn  VIIL=HFF«MFN.ilH*440,KJC;L  I  B.m  SFs(  NEK.*eeF  I  , 

•  //  OsN*MFN,llHAA90.KJC.FLnT  .IHITPIIT  ,SMACF»  tCYL,  (1,1)  ,RLSF|, 

//  (ICH«(RFcFM«FH,LRFCL»BO.im.KSF7FsAlgO  )  _ 

•  //nATA.  iNi»iiT  no  * 
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Program  VI 

*Program  to  compute  relative  HF  power  above  isotropic  on  a  70km 
altitude  plane.  Frenquency=3. 1 7  MHz 

*note:  Function  subprograms  AF  (Theta,  Phi)  and  ELF  (Theta,  Phi) 
are  not  shown.  They  are  the  same  as  in  program  IV. 


mu m j*  ij ii  •  .’Di •(» ‘••ii 
u  i’ )  jA  1  i  JO  -  I  »A  •  I ,  j  j » j  j  *  n » s  .:s 
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Program  VII 

*Program  to  compute  relative  HF  power  above  isotropic  on  a  70km 
altitude  plane.  Frequency=5.lMHz. 

♦note:  Function  subprograms  AF  (Theta, Phi)  and  ELF  (Theta,  Phi) 
are  not  shown.  They  are  the  same  as  those  in  program  V. 


Pt h; '<«. 'll <>  icie b/ '(Sit ^ / •j'iM ^ I (ibii/'js'/t <f I o;.Bi ^ I uhHi"*,*/! .. lui.b/ i(.  -hi 'IS*-* ('ll (.M/  I u'irfi  »sbi  c: Ipi  lu ■ 


//rc^  f.  i-isr-  l<'L».r-.-Tri  '>  v.-i'vn) 

,  ,  L;t  =  '  i'i“  1. 1  .u  .  PI  rT.  f  "1  nn 

//  IXF''  ‘i-cr 

/•.II  E!  •  *ir'  =  riFEP 


ii j.i *.1  ih j  1  - < n  '.if 


Of.i  /  f  I  Jbut ‘#sn{2  I  2  ICtiyi  4S»il  ?:  lOtbf^SM  t  llroUi#S**|  I’Uf  .li  !  1 1  t k  Jl  C'l.M  t'^SM  e  1 1  i.rt  iS* 


l,  6  f,  I 


Program  IX 


Program  to  compute  relative  magnetic  field  strength  of  an 
observation  point  for  an  ionspheric  ELF/VLF  current  element 
5.1  MHz  HF  pattern. 


)'IM.  J'lS'.lM-ltllJ  JUSi.  71111 IJ  1'l'if)7«1  5  1  2  IllllTHin  1  2  lll'.TH'in  V  IIS'  7  810 


:  I"  i.i«d  li  !..  -  u 

4  i  Mvi  M  •  d  1  it  J  i  =  O 
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APPENDIX  II 

USE  OF  CHEBYSHEV'S  POLYNOMIALS  TO  SIMPLIFY  ANTENNA  FACTORS 

The  Chebyshev's  polynomials  are  the  solution  to  the  Chebyshev 
differential  equation  (A-1).  The  solution  has  the  fonn  of  equation 
(A-2) 

(1  -  x2)^— ^  +  n^y  »  0  (A-1) 

dx2  dx 

With  a  recursive  formula  given  in  equation  (A-3)^^^^ 

Tnj(x)  =  cos  (m  coa“^x)  m  >  0,  1x1  <  1  (A-2) 

T„+i(x)  -  2x  T„(x)-T„_i(x)  (A-3) 

Equation  (A-4)  and  (A-5 )  follow  from  substitution  of  0  and  1  for 
m  in  (A-2) 

To  »  1  (A-4) 

Tj  ■  X  (A-5) 

By  using  the  recursion  relationship  (A-3),  the  following 

polynomials  are  obtained: 

T2(x)  *  2x2  _  1 
T3(x)  *  4x3  _  3x 
T4(x)  -  8xA  -8x2+1 

T5(x)  *  16x3  -  20x3  +  5,^ 

T6(x)  -  32x6  _  48x4  +  i8x2  -  1 
T7(x)  "  64x2  _  112x3  +  56x3  _ 

T8(x)  •  128x8  -  256x6  +  isox^  -  32x2  +  1 

Let  "w"  equal  "co3~^x”  then  "x"  is  equal  to  "cos  w. " 


Substitution  for  "x"  in  equation  (A-2)  gives  (A-6). 
Tn(co8  w)  *  cos  (m  wl 


(A-6) 
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By  using  equation  (A-6)  with  the  polynomials,  trigonometric 
identities  can  be  found  for  expressing  "cos  (m  w)"  or  "sin  (m  w)"  in 
terms  of  "sin  w"  and  "cos  w."  For  example,  using  "T2(x)"  and  letting 
"x"  equal  "cos  w" ,  equation  (A-7)  is  obtained. 

T2(cos  w)  =  cos  (2w)  =  2  (cos^  w)  -  1  (A-7) 

.Expressions  for  sin(m  w)  can  be  obtained  by  taking  the  derivative 
of  the  "cos  (m  w)"  identity.  Equation  (A-8)  was  obtained  by  taking 
the  derivative  of  (A-7). 

sin  (2w)  =  2  (cos  w)  sin  w  (A-8) 

To  obtain  the  antenna  factors  in  the  form  of  equation  (1-13)  and 
(1-14),  equation  (1-7)  must  be  expanded.  Equation  (A-9)  is  the 
expansion  for  the  8-element  case.  Let  w  =/3(d/2)  sintf  . 

8-1 

AF  *  ^  cos  (m/3d/2  sin  0)  *  cos(w)  + 

in»l  ,3,5. . . 

cos  (3w)  +  cos  (5w)  cos  (7w)  (A-9) 

Use  the  Chebyshev's  polynomials  to  find  trigonometric  identities  to 
reduce  (A-9)  into  an  equation  in  terms  of  powers  of  cos  w.  These 
identities  are  given  in  equation  (A-10). 

Tj  =■  cos  w  *  cos  w 

T3  »  cos  3w  *  4cos^  w  -  3  cos  w  (A-10) 

T5  ■  cos  5w  =•  Ibcos^w  -  20cos^w  +  5cos  w 

T7  =•  cos  7w  =•  64cos^w  -  112cos5w  56cos^w  -  7cos  w 

Substituting  (A-10)  into  (A-9),  the  expression  for  the  antenna 

factor  becomes  (A-11). 

AF  *  64cos^w  -  96co85w  +  40cos -w  -  4cosw  (A-11) 
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An  identity  for  "cos(8w)"  can  be  found  from  polynomial  "T3."  By 
taking  the  derivative  of  "cos  (8w)",  an  identity  for  "sin  (Sw)"  can  be 
found,  (A-12). 


sin  8w  =  1024  cos^w  sin  w  -  1536  cos^w  sin  w  +  640  cos^w  sin  w  - 
64  cos  w  sin  w  =  16  sin  w  (64cos^w  -  96  cos^w  +  40  cos^w  -  4cos  w) 


(A-12) 


Equation  (A-12)  can  be  rearranged  into  (A-13). 


Isin8w_,,  7  7 

—  — : -  =  64  cos'w  -  96  cos-'w  +  40  cos-^w  -  4  cos  w 

16  sin  w 

(A-13) 

Equation  (A-14)  can  be  obtained  by  subtituting  equation  (A-13) 
into  (A-11). 

1  3  in  8w 


AF 


(A-14) 


16  sin  w 

This  is  identical  with  equation  (1-14)  for  the  8-element  array 
antenna  factor  with  the  exception  of  the  1/16  constant  in  (A-14).  The 
constant  can  be  neglected  at  this  point,  due  to  the  fact  that  when  the 
directive  gain  is  calculated,  the  AF  will  be  normalized. 

A  similiar  calculation  can  be  performed  to  obtain  the 


4-element  array  factor,  (1-13). 
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APPENDIX  III 

REVISED  ARECIBO  HF  ANTENNA  ARRAY  GEOMETRY 


After  completion  of  the  work  described  in  the  main  body  of 
this  report,  a  preliminary  copy  was  reviewed  by  the  Arecibo 
Observatory,  At  this  time  the  Arecibo  Observatory  made  available 
additional  information  about  the  HF  array.  The  following 
differences  were  reported  between  the  model  used  in  the  main  text 
and  the  actual  A.O.  array: 

1.  The  T  for  the  antenna  is  .774. 

2.  The  pyramid  structure  is  elevated  five  feet  above 
ground. 

3.  The  feed  point  of  each  face  is  capacitively  loaded. 
This  appendix  presents  the  information  provided  by  the  A.O.  and 
discusses  its  possible  effects  on  the  work  presented  in  the  main 
body  of  this  report. 

The  T  of  the  NPLA  is  .774.  It  is  obtained  by  taking  the 
ratio  of  the  lengths  of  two  consecutive  elements  of  a  face  on  the 
same  side  of  the  feed  line.  Table  III-l  provides  a  list  of  the 
element  lengths  and  the  length  of  the  feed  line  to  that  element. 
The  even  numbered  elements  represent  the  dimensions  for  the 
scaled  faces  of  the  pyramid  structure. 

Figure  (III-l)  provides  a  view  of  the  structure  at  the 
vertex  of  the  pyramid.  It  was  concluded  from  this  figure  that 
height  of  the  vertex  of  the  pyramid  was  1.524  m  (60  inches). 
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Element 

Element  Length 

Feed  Line  Length 

Number 

(m) 

(m) 

1 

25.000 

35.704 

2 

23.454 

33.49 

3 

22.001 

31.42 

4 

20.638 

29.47 

5 

19.361 

27.65 

6 

18.163 

25.94 

7 

17.035 

24.33 

8 

15.984 

22.82 

9 

14.993 

21.41 

10 

14.064 

20.01 

11 

13.195 

18,84 

12 

12.378 

17.68 

13 

11.610 

16.58 

14 

10.891 

15.55 

15 

10.217 

14.59 

16 

9.586 

13.69 

17 

8.992 

12.84 

18 

8.434 

12.05 

19 

7.913 

11.30 

20 

7.422 

10.60 

21 

6.902 

9.94 

22 

6.532 

9.33 

23 

6.126 

8.75 

24 

5.749 

8.21 

25 

5.392 

7.70 

26 

5.060 

7.22 

27 

4.746 

6.78 

28 

4.450 

6.36 

29 

4.176 

5.96 

30 

3.917 

5.59 

III-l.  NPLA 

Element  and  Feed  Line 

Lengths  Provided  by 

Even 

element  numbers  are  for  the  scaled  faces. 

T  -  . 774. 
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Figure  III-l  Vertex  of  pyramid  for  Arecibo  Observatory  HF  non-planar  log-periodic  array 


Figure  (III-2)  shows  a  sketch  of  a  feed  point  of  one  of  the 
faces  of  the  pyramid.  The  wires  detailed  by  "A"  are  in  addition 
to  previously  used  geometry.  All  of  the  conducting  elements  of 
the  antenna  consist  of  three  #12  twisted  steel  wires  with  10% 
aluminum  coating. 

The  additional  information  on  the  array  elements  was 
combined  into  a  new  geometry  description  for  AMP.  Discrepancies 
exist  within  the  new  information  provided.  Figure  (III-l) 
depicts  the  height  of  the  vertex  at  five  feet,  while  figure 
(III-2)  shows  a  six-foot  height.  Table  III-l  and  figure  (III-2) 
give  different  distances  between  the  feed  point  and  the  bottom 
element.  It  is  unclear  whether  figure  (III-2)  represents  a 
scaled  face  or  unsealed  face  of  the  pyramid  and  whether  the 
dimensions  of  the  additional  wires  are  also  scaled  between  the 
two  sets  of  faces.  For  the  new  geometry  description  to  be  used 
with  AMP  ...le  following  geometry  was  decided  upon: 

1.  The  height  of  the  vertex  is  1.524  m. 

2.  The  lengths  of  all  the  feed  lines  to  the  elements 
are  all  taken  from  table  III-l  (including  the  length 
to  the  bottom  element). 

3.  Figure  (III-2)  was  taken  to  be  an  unsealed  face  of  a 
pyramid.  The  144"  dimension  was  taken  to  be 
correct  and  the  other  dimensions  were  adjusted  to 
correspond  to  table  III-l. 

4.  All  wires  and  dimensions  from  the  unsealed  face  were 
scaled  by  the  fourth  root  of  t  (t  *  .774)  for  the 
scaled  face.  This  includes  the  additional  wires 
shown  in  figure  (III-2). 

The  final  AMP  geometry  deck  incorporating  these  changes  is 
given  in  figure  (III-3).  The  computer  results  of  the  power  gain 
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Figure  III-2  Capacitvely  loaded  feed  region  for  one  face  of  pyramid  element  in 
Arecibo  Observatory  HF  heating  non-planar  log-periodic  array 
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£or  the  two  cases  of  3.17  MHz  and  5.1  MHz  are  shown  in  figures 
(III-4)  and  (III-5).  In  these  figures  the  "x"  denotes  the 
results  of  the  new  modified  geometry  given  in  this  appendix,  and 
the  denotes  the  results  of  the  "old"  geometry  used  in  the 
main  body  of  this  report.  Figure  (III-4)  shows  the  power  gain  as 
a  function  of  "phi"  for  selected  constant  values  of  "theta." 

These  are  the  same  values  of  "theta"  used  in  figure  (1-6). 

Figure  (III-5)  is  a  plot  of  power  gain  as  a  function  of  "theta" 
for  selected  values  of  "phi."  The  selected  values  of  "phi"  are 
the  same  values  of  "phi"  used  in  figure  (1-7)  plus  additional 
values  corresponding  to  the  "x"  and  "y"  axis  (i.e.,  (p  =  0",  90°, 
180°,  270°). 

Examination  of  figures  (III-4)  and  (III-5)  leads  to  the 
conclusion  that  the  changes  in  geometries  between  the  two  cases 
result  in  only  a  small  difference  in  power  gain  for  small  values 
of  "theta"  and  large  differences  for  large  values  of  "theta." 

It  is  necessary  to  determine  what  effect  the  new  geometry 
will  have  on  the  results  of  the  main  body  of  this  report.  The 
elemental  power  gain  for  0  <  50°  is  approximately  equal  for  the 
two  geometries.  Since  a  large  portion  of  the  radiated  power  is 
contained  in  this  region,  it  is  expected  that  the  directivity  of 
the  total  array  would  remain  approximately  the  same.  The  results 
shown  in  figures  (1-10)  and  (1-11)  should  be  approximately  the 
same  for  0  <  50°  but  significant  differences  could  occur  for 


Q  >  50°. 


175 


Figure  III-4a 


Power  gain  vs.  phi  for  constant  theta.  Comparison 
of  "new"  and  "old"  heating  array  element  geometrv. 
Frequency=  3.17  MHz 


Figure  III-4b  Power  gain  vs.  phi  for  constant  theta.  Comparison 
heating  array  element  geometry.  Frequency*  5.1  MH 
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'old'  heating  array  element  geometry.  Phi  =  90'^ 
uency  =  3.17  MHz 
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Figure  III-5a.5  Power  gain  vs.  ’theta’  for  constant  ’phi’.  Comparison  of  'new 
and  ’old’  heating  array  element  geometry.  Phi  =  40° 

Frequency  =  1.17  MHz 
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Figure  III-5b.4  Power  gain  vs.  'theta*  for  constant  'phi'.  Comparison  of  'new 
and  'old'  heating  array  element  geometry.  Phi  =  270° 

Frequency  =  5.1  MHz 
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Figure  IH-5b.6  Power  gain  vs.  'theta'  for  constant  'phi'.  Comparison  of  'new 
and  'old'  heating  array  element  geometry.  Phi  =  250° 

Frequency  =  5.1  MHz 
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Of  particular  interest  is  the  differences  in  directive  gain 
along  the  "x"  and  "y"  axis.  The  ELF/VLF  source  array  orginates 
from  the  HF  antenna  response  along  these  lines.  Table  III-2 
shows  the  value  of  "theta"  to  the  center  and  furthest  edge  of  the 
ELF/VLF  source  region  (shown  in  figures  (1-15)  and  (1-16))  which 
is  furthest  from  the  origin. 


Frequency 

Source 

X 

Location 

y 

Length 

Width 

Center 

Edge 

3.17  MHz 

63  km 

0 

24  km 

42“ 

47“ 

3.17  MHz 

0 

38  km 

— 

7  km 

29“ 

31“ 

5.1  Mhz 

60  km 

0 

54  km 

— 

41* 

51“ 

5.1  MHz 

0 

65  km 

28  km 

43“ 

49“ 

Table  I I 1-2. 

Value 

for  0  to  the 

Source  Regions  Furthest 

from  Origin. 


Since  all  the  ELF/VLF  sources  are  located  in  a  region  where 
"theta"  is  less  than  SO”,  it  is  concluded  that  the  new  geometry 
will  not  significantly  affect  the  zero  approximation  ELF/VLF 
array  model.  The  plots  of  relative  field  intensity  versus 
ELF/VLF  frequency  should  remain  essentially  the  same.  The 
fundamental  conclusion  of  the  main  body  of  the  report  that  the 
ELF/VLF  frequency  response  is  affected  by  the  geometry  of  the  HF 
heating  antenna  pattern  remains  intact. 
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